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PREFACE 



In 1942, the author published a monograph entitled 
"Heat Transfer by Infrared Radiation in the Atmos- 
phere," which was rather favorably received by the 
meteorological public. The war terminated this par- 
ticular work, and during it and thereafter the author 
became engaged in other scientific pursuits not directly 
related to atmospheric radiation. In the meantime, a 
vast literature has grown up around the subject, dealing 
with laboratory experiments as well as theoretical 
analysis, not to speak of the applications to the 
atmosphere itself. From the fall of 1950 until 1954, 
the author, while at the University of Utah, was 
engaged in a project sponsored by the Air Force 
Cambridge Research Center whose aim it was in 
part to extend the methods used earlier for water vapor 
so that they were applicable to carbon dioxide and 
ozone. Very extensive new laboratory data on all 
three optically active atmospheric gases have become 
known in the last seven or eight years. This finally 
suggested a new edition of the earlier monograph, 
and the task was begun in the summer of 1957 at the 
Scripps Institution of Oceanography, again under the 
sponsorship of the Air Force Cambridge Research 
Center. 

The present monograph is a second edition of the 
earlier one in the limited sense that the methods used 
to determine radiative fluxes in the atmosphere have 
been retained with only very minor changes. Approxi- 
mations much more refined but also much more 
cumbersome to handle than the linear techniques of 
the 1942 work have been omitted. Instead, radiation 
tables have been calculated on the basis of the avail- 
able data, practically all of them fairly new, for the 
15-/z C0 2 band, the combined 9.6-/* and 9.0-/x ozone 
bands, the 6.3-ju water-vapor band, the rotational 
water- vapor band, and finally for the overall water- 
vapor spectrum. The only difference, in principle, 
from the earlier treatment lies in the use of a linear 
pressure reduction throughout, in place of the former 
square-root reduction. In the text, the review of the 
basic theory of band structure and band absorption 
has been somewhat enlarged. On the other hand, 
Part III of the earlier monograph which dealt mainly 
with radiation measurements in the atmosphere itself 
has been omitted, partly because the new tables are 
based on laboratory data, but largely because the 



literature on the subject has swollen to such dimen- 
sions that it would be impossible to do justice to it 
within the confines of this monograph. 

Owing to many other obligations, the author has 
been unable to keep up in detail with the proliferating 
literature on atmospheric radiation. A purely formal 
bibliography based on a search of the abstracting 
literature seemed out of place in a work of this limited 
character. We have therefore restricted the bibliog- 
raphy to papers that are actually quoted in the text. 
It is possible, though not too probable, that an occa- 
sional paper which has a direct bearing on our work 
might have been overlooked, and the author apologizes 
in advance for such omissions. 

Another restriction of the present monograph con- 
cerns the fact that we give only tables and no 
nomograms (charts). This is based on the view that 
a work like the present by its very nature addresses 
itself not so much to the practical forecaster as to the 
scientific investigator of the atmosphere. The latter 
will almost invariably spend some time and effort on 
the subject of radiative transfer if he becomes engaged 
in it at all. Given the very numerous graphical or 
numerical schemes that may be used for this purpose, 
depending on individual tastes, we felt it to be prefer- 
able to leave the construction of nomographs to the 
investigators concerned. With modern methods of 
reproduction, this is a relatively simple and cheap 
task that can be performed quite rapidly. 

The very extensive computational work, including 
the involved two-dimensional smoothing of the radia- 
tion tables, was done by Miss Margaret Culbertson, 
M.A. Without her consistent accurate performance, 
this monograph could not have been written. 

We are greatly indebted to the Thermal Radiation 
Laboratory of the Air Force Cambridge Research 
Center for their financial support, and in particular 
to Drs. John N. Howard, J. I. King, and P. R. Gast 
for their personal encouragement. This assistance and 
encouragement has been invaluable in completing 
the work. 

Walter M. Elsasser 

Scripps Institution of Oceanography 
La Jolla, California 
February 1960 
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I. INTRODUCTION 



This monograph contains numerical tools required 
to compute approximately the effects of radiative 
transfer in the atmosphere for the far infrared region 
of the spectrum. We shall not deal with the transmis- 
sion and scattering of solar radiation, but only with 
the long-wave radiation that is emitted, and of course 
also absorbed, by the earth's surface, by clouds, and 
by the atmosphere itself. This part of the spectrum 
extends from a wavelength of about 4 microns to near 
the microwave range at a fraction of a millimeter 
wavelength. The main agents of transfer are two 
water-vapor bands and one band of carbon dioxide. 
In addition, there are some bands of carbon dioxide 
of minor importance and a band or rather a cluster of 
bands of ozone. 

It is well established that the net effect of infrared 
radiation transfer is almost invariably a cooling. 
Only in a few rather exceptional circumstances does 
there occur a small amount of heating. Radiative 
cooling takes place at the ground, especially in clear 
weather, at cloud tops, and in the free atmosphere. 
The average radiative cooling of the free atmosphere 
is not large; it is of the order of one to two degrees 
per day, but, since this goes on all the time, the cumu- 
lative effects are of major importance for the dynamics 
of the atmosphere. The present monograph provides 
data from which all these effects of radiative cooling 
can be calculated simply and with a reasonable de- 
gree of accuracy at any given location, provided the 
temperature and moisture distribution of the at- 
mosphere are known from a sounding. 

Far infrared spectroscopy is now about half a 
century old and highly developed. The principal 
bands, especially those of simpler compounds such 
as H 2 0 and C0 2 are exceedingly well investigated; 
almost all of their lines have been measured and have 
been interpreted in quantum-mechanical terms. More- 
over, even the lesser bands are now rather thoroughly 
known and their molecular mechanism understood. 
Unfortunately, this does not solve the problem of 

1 Present affiliation : Chairman, Department of Physics, Uni- 
versity of New Mexico, Albuquerque, New Mexico. 



atmospheric transfer. Each band consists of hundreds 
of lines, and the total transfer is obtained by combining 
the contributions from all of these lines. Proposals 
have frequently been made to deal with this problem 
by means of large-scale computing machines. This can 
no doubt be done, but it would involve the manipula- 
tion of a vast amount of primary data on the part of 
the investigator, and it would even then be difficult 
to know whether all the significant data and effects 
have been taken into account. While electronic com- 
puting machines have so far seldom been used for this 
purpose, Kaplan (1950, 1952) has carried out elaborate 
manual computations for the 15-ju carbon-dioxide 
bands. The results are in good agreement with observa- 
tions; but even after smoothing over the individual 
lines there remains a good deal of cumbersome detail, 
which makes it difficult to generalize these results to 
as broad a variety of conditions as we meet in the 
atmosphere. 

It is hard to see how a rigorous solution of the trans- 
fer problem can be found, applicable to a large variety 
of atmospheric conditions, given the complexity of 
the spectrum and the intricate way in which the line 
intensities depend on temperature and the line shapes 
on pressure. On the other hand, there exist approxima- 
tions which in the majority of cases are satisfactory. 
The nature of these approximations is described in 
detail in section II; their physical background is 
described in section III. They follow in a straight- 
forward way from the basic formalism of the radiative- 
transfer equations together with typical structural 
characteristics of molecular bands. Some important 
points regarding these approximations should be 
brought out here. 

In the first place, averaging over spectral intervals 
containing several or many lines is carried out mathe- 
matically at an early stage of the formalism, by means 
of rather straight-forward theoretical models. In 
terms of these averages, solutions of the transfer 
equations can be found that are rigorous for a schemat- 
ical band and should be very nearly correct for real 
bands. The only features of the bands entering there- 
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after are smoothed-out band contours whose trans- 
mission characteristics can be expressed in terms of a 
very few numerical parameters. These parameters 
lend themselves to direct observational determination 
in a manner to be indicated. 

Band-spectroscopic observations can be carried out 
in two ways. The traditional ambition of the spectros- 
copist has been to work with as narrow a slit as 
possible in order to resolve all the spectral lines, even 
the very weak ones. This is, however, not so advan- 
tageous when dealing with problems of atmospheric- 
radiation transfer. Here, we want a slit wide enough 
so that we do not see the individual lines at all but, 
instead, measure directly a smoothed-out band con- 
tour. For a long time, no such measurements were 
available, but in the last seven or eight years system- 
atic investigations of far infrared bands with wide- 
slit techniques have been undertaken, using specially 
designed absorption cells to simulate or at least 
approximate the long optical paths encountered in the 
atmosphere. For the main bands of water vapor and 
carbon dioxide, this type of work was done at the Lab- 
oratory of Astrophysics and Physical Meteorology of 
Johns Hopkins University and at the Physics Depart- 
ment of Ohio State University. Many of the measure- 
ments were made with somewhat differing techniques at 
both places, providing a welcome mutual check. 
Furthermore, these measurements were made in the 
presence of various quantities of air in the absorption 
cell, allowing one to study directly the very important 
effects of extraneous pressure upon the absorption 
characteristics of the far infrared bands. This labora- 
tory work is described in section IV. It is not possible, 
unfortunately, to cool the extremely bulky absorption 
cells below room temperature, and so the extrapolation 
of the band absorption to the low temperatures of the 
upper atmosphere has to be done by theoretical means. 
There, again, we have used an approximative averag- 
ing method and did not go back to the individual lines. 

The tables given here are exclusively based on 
laboratory measurements of the above-mentioned 
type that simulate atmospheric conditions. On the 
one hand, direct reference to finer structural details 
of the bands is avoided ; on the other hand, there is also 
no reference to data gained in the atmosphere itself. 
Earlier experience with water-vapor transmission indi- 
cates that the actual radiative fluxes will be closely 
reproduced, probably to within a few per cent in most 
cases, by these tables ; nevertheless, the ultimate check 
will have to consist of extended measurements under 
a great variety of atmospheric conditions, to be com- 



pared with the predictions obtained from the tables. 

The tables form the basis of radiative-cooling cal- 
culations. They may also be used to compute the total 
radiation leaving the earth which is of interest on 
general geophysical grounds and to determine the 
flux received by a satellite. All the numerical work 
may be done once and for all except for one final 
numerical integration which involves the temperature 
distribution of the atmosphere as well as the distribu- 
tion in height of the optically active gases, and there- 
fore it varies with conditions, both in space and in 
time. Very early, Mugge and Moller (1932) developed 
a nomogram on which this integration could be done 
graphically. It was further developed by Elsasser 
(1942) and has become known as a 4 'radiation chart." 
We retain this term, although we present the data 
only in tabular form. Many workers prefer to do the 
integrations numerically rather than graphically. Also, 
we find that, if a nomogram is desired, it can be inked 
by a draftsman and reproduced by some conventional 
process in very short order and with a minimum of 
expense. The term ' 'radiation chart" appearing in the 
following pages designates a table with two entries, 
temperature, T, and optical thickness in the vertical, 
w, so that, when the empirical relationship between 
the two, say u(T) y is given for any particular atmos- 
phere, the radiative flux may be obtained by integra- 
tion along the path defined by u(T). This technique 
is described in section V. 

The question has often been asked whether radia- 
tive-cooling calculations could be used for purposes 
of forecasting. We believe that the calculations are, 
as a rule, prohibitively difficult; moreover, radiative 
cooling is not usually strong enough to affect the 
local dynamics of the weather over a relatively short 
period. Apart from such local effects as nocturnal 
cooling of the ground, radiative cooling is a "secular" 
effect in the dynamics of the weather. Its variation 
with altitude, latitude, season, etc. should be of great 
meteorological importance; these effects can be in- 
vestigated systematically by means of the data pre- 
sented here. So far as altitude is concerned, our meth- 
ods should be reasonably accurate up to a height of, 
say, 25 to 30 km but should give at least correct 
orders of magnitude to, say, 40 km. This is discussed 
in section 3.2. A number of investigators have carried 
out computations of radiative cooling based on earlier 
and, we believe, already fairly accurate data. The 
material contained in this monograph should permit 
a further check of such work and an extension to a 
greater variety of atmospheric conditions. 
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II. PRINCIPLES OF RADIATIVE TRANSFER 



1. Kirchhoff's law, Planck's law, and the equation of 
transfer 

The problem encountered in determining the infra- 
red fluxes in the earth's atmosphere is much simpler 
than the general radiative transfer problem occurring 
in geophysics or astrophysics. There are two essential 
simplifying facts which hold to a very high degree of 
approximation. The first is the absence of scattering. 
A parcel of air containing optically active gases will 
partly absorb and partly transmit incident radiation, 
but it will not appreciably scatter the transmitted 
radiation. In the second place, every particle of air is 
effectively in local thermodynamical equilibrium. 
That is, the excitation of molecular levels and the 
radiation ensuing therefrom are the same as they 
would be if the parcel were in an enclosure in thermo- 
dynamical equilibrium at a given temperature, T. 
The condition of the parcel, including its emitted 
radiation, can therefore be described by density, 
temperature, pressure, and the relative concentrations 
of the optically active substances in the parcel; the 
radiation issuing from the parcel is, in particular, 
independent of the characteristics of the radiation 
incident upon the parcel. In the next chapter, we shall 
briefly return to these conditions and their physical 
background ; for the time being, we shall admit them 
as valid. 

A surface which absorbs all radiation of any given 
wavelength falling upon it is called completely opaque. 
A surface which is completely opaque at all wave- 
lengths is called black. An enclosure surrounded by 
black walls, the whole kept at a constant temperature, 
T 9 is filled throughout by black-body radiation of 
temperature T. This is true no matter what else is 
inside the enclosure. If local thermodynamical equi- 
librium prevails in the enclosure, Kirchhoff's law can 
readily be proved. Consider a piece of matter inside 
the enclosure, and consider a given direction of travel 
for the radiation. Incident upon the piece, there is in 
the given direction black-body radiation of intensity, 
say, Jfc(X) for any given wave length, X. The fraction 
A of this radiation may be absorbed. Let E be the 
radiation emitted by the piece. Then, in equilibrium, 

h = (1 -A)I h + E, 

expressing the fact that the total radiation leaving the 
piece is again black. Therefore, 

E/A = / 6 , (1) 

which is Kirchhoff's law. This is valid at any specified 
wavelength. 



In this paper, we shall use as a unit the wave 
number, which is the reciprocal of the wavelength; 
we designate it by v throughout, and 

„ = l/x = v'/c, (2) 

where v' is the conventional frequency and c the 
velocity of light. 

We shall now indicate Planck's law of black-body 
radiation together with the numerical constants that 
have been used in the present tables. We encounter 
two situations of interest. On the one hand, we con- 
sider unidirectional beams of radiation of intensity 
h(v) per unit area as emitted by a black body; on 
the other, we are interested in the black-body flux B v 
that emerges from a unit area of a black surface into 
all directions of the hemisphere. The two are con- 
nected by the simple relationship 

B v = wh(v) (3) 

which will be proven later on. It is customary to ex- 
press Planck's law in terms of fluxes. Then 

pv z 

B > = ~r (4) 

where, in terms of fundamental constants, 

p = 2whc 2 , q = hc/k 

(h, Planck's constant, k Boltzmann's constant) and 
numerically (figures from American Institute of 
Physics Handbook, McGraw-Hill, N. Y., 1957; see, 
also, Cohen et al., 1955) 

p = 3.7413 -lO" 5 erg/cm 2 sec 

q = 1.4389 cm deg. (5) 

The absolute temperature is related to the centigrade 
scale by 

T = t + 273.16. 

We shall frequently use the derivative with respect 
to the temperature of the black-body function. By 
letting 

x = qv/T, 

this can be written 

dB v pT 2 x 4 e x 

where 

p/q* = 1.25585 • 10" 5 , 
using the values from (5). In our tables, we have 
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Fig. 1 . Shape of spectral line. 

adopted as a unit one gcal per sq cm per day. In 
these units, the factor becomes 



general a very complicated one. The far infrared bands 
which are responsible for radiative transfer in the 
atmosphere appear, on sufficient spectral resolution, 
to consist of very numerous spectral lines, and we 
shall proceed by first considering the absorption of a 
single spectral line. We shall then study band patterns 
made up of many individual lines. 

The shape of spectral lines can depend on a variety 
of factors, but in most of the atmosphere the far infra- 
red lines have the so-called Lorentz shape. We shall 
designate the integrated intensity of a line by 5. 
Then the distribution of absorption over the line 
is expressed by the formula (fig. 1) 



p/q* = 2.5915 -10- 8 (gcal/cm 2 -day). 



(7) 



For the sake of completeness, we mention the figure 
for the total black-body radiation (Stefan-Boltzmann 
formula) 



B = f B4v, 
Jo 

B = 1.1696- 10- 7 r 4 (gcal/cm 2 <day) 
dB 

= 4.6784- IQrT* 

dT 



(8) 
(9) 



In order to derive the equation of radiative transfer, 
we now return to the case of a straight-line beam, 
leaving the integration over the hemisphere for later. 
We shall use u to designate optical thickness (mass 
of optically active gas per sq cm). Let k be the 
absorption coefficient, defined by the relation 



dl a = - kldu ; 



(10) 



that is, in a layer of optical thickness du, the fraction 
kdu of the incident intensity is absorbed. By Kir- 
chhoff's law (1) the amount emitted by the layer is 



and hence 



die = klbdu, 
dl = - k(I - h)du, 



(ID 
(12) 



S a 
k(v) = V2 , 9 - 



This is normalized, so that 



kdv — S. 



(14) 



(15) 



The parameter a is called the half-width of the line ; 
it is a function of pressure and, to a lesser degree, of 
the temperature. For a constant temperature, a is 
directly proportional to the total pressure, which in the 
atmosphere is essentially the air pressure, since the 
partial pressures of the optically active gases are quite 
small. 

We now define the transmissivity of a spectral in- 
terval Av as 

T = _JL f Uv =— f e-*™*dv, (16) 
I^Av J Av J A „ 

by (13). By substituting here the expression (14), 
we get 

Sail/ it 



If f Sau/ir I 



(17) 



which is the equation of radiative transfer for the one- 
dimensional case in the absence of scattering. 

2. A single spectral line 

Next we shall be concerned with straight absorp- 
tion, expressed by the integral of (10)— 

I = he- k \ (13) 
Here, k is a function of v, and in band spectra it is in 



If we take the interval wide enough, we may use - =o 
and +oo as limits of integration. We introduce the 
abbreviation 

x = Su/2ira. (18) 
By substituting the variable, 

5 v — vq 

tan - = , 

2 a 



(17) becomes 

Av J- 



p—x cos s—x 



d ( S \ a 

— I tan - I ds. 
ds V 2/ 



By an integration by parts and some further manipu- 
lations, this finally assumes the form 

A = 1 - r = (2ira/Av)xe-*ZJo(ix) - i/i (**)]. (19) 



I 
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Here, A will be called the absorptivity, Jo and Ji are 
Bessel functions of order zero and one, respectively. 
This formula is due to Ladenburg and Reiche (1911). 

There are two important limiting cases. When x is 
small, the power series for the Bessel functions yields 

AAv — 2irax = Su. (20) 

This is the region of linear absorption, the absorption 
being proportional both to line strength and to the 
optical path. The condition of validity of this ap- 
proximation is clearly x « 1 which, by (18) and (14), 
means 

k(va)u <K 1 ; 

i.e., for this approximation to hold, the absorptivity 
must be small at the line center; it is then a fortiori 
small in the wings of the lines (fig. 2). 

In the opposite extreme, x is large. The asymptotic 
development of the Bessel functions then gives 

AAv = 2V5W (21) 




Fig. 3. Periodic band pattern. 



The absorption is now proportional to the square root 
of the line intensity and of the optical thickness. 
One can readily verify that, as fig. 2 shows, the central 
part of the line has been completely absorbed away, 
and further absorption can occur only in the wings of 
the line. 

Formula (21) may be obtained more directly by 
neglecting the term a 2 in the denominator of (14). 
If this simplified form for k is substituted into (16), 
one may obtain (21) by an elementary integration. 

For a band containing many lines, (20) leads to 

AAv — £ SiU, 

the summation extended over all the lines in the in- 
terval. This is a good approximation when each term 
of the sum is sufficiently small, since mutual over- 
lapping of neighboring lines is unimportant for weak 
absorption. On the other hand, for the strong absorp- 
tion case expressed by (21), overlapping of neighboring 
lines in a band becomes of paramount importance, 
and the total absorption will always be less than the 
sum of the contributions (21) of the individual lines. 
Clearly, A cannot exceed unity, whereas the preceding 
expressions for A are unbounded and therefore in- 
applicable to sufficiently thick layers. The next de- 
velopments are concerned with taking into account 
the mutual overlapping of lines in a band. 



3. Periodic band pattern 

A simple and in some cases fairly realistic model of a 
band suggests itself readily: we allow a single line of 



the form (14) to repeat itself periodically (fig. 3). 
Thus, we let 



n=H- 00 



Sa/ic 



(v - ndY + a 2 



(22) 



where d is clearly the distance between adjacent lines. 
It is possible to express (22) in terms of an analytical 
function owing to the fact that, if such a function 
has only single poles, it is uniquely defined by these 
poles (by Mittag-LefHer's theorem ; see Whittaker and 
Watson, 1940, Ch. 7). We find that (22) is equiva- 
lent to 



S 



sinh 0 



d cosh ft — cos s 



(23) 



where 



$ = 2ra/d, s = 2wv/d. (24) 
The transmissivity (16) may now be written 



1 r+* 
T = — I e-*w*ds. 

2 IT J—* 



To evaluate this, form 

dr 
du 



= f ke~ ku ds. 

2tJ 



We introduce the variable 

cosh cosh s — 1 

cos <p = — — 

cosh j3 — cosh s 
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and find 
dr S ^ 
du 



s r 2ir r Sm ~\ 

= - — I dip exp — . (coshff+cos ^) 
2?r^Jo L d-sinh/3 J 



with 



d • sinh /? 
= - (S/d)e~* coth ^/o(^/sinh 0), 

3; = 



(25) 



By integrating the last expression with suitable 
boundaries, 



= f e -v ooth 0Jo(iy/smhp)dy. 

Jy 



(26) 



We shall consider limiting cases of this formula. The 
parameter measures the ratio of line width to the 
distance of neighboring lines, by (24). For large 0 
(that is, strongly overlapping lines), we have coth 0= 1 
and Jo = 1 ; hence, 

r = e~* t 

and the band exhibits the behavior of a continuous 
absorption spectrum. This case, however, is not 
realistic for the atmosphere. We should note that the 
degree of overlapping does not only depend on /3 but 
also on the thickness u of the layer. The case of large /5 
corresponds to a "smearing out" of the absorption 
coefficient even if u is so small that the absorption of 
the individual line is represented by the "linear" 
formula (20). On the other hand, overlapping may 
occur in thick layers (large u) where the individual 
lines absorb only in the wings. This is the normal case 
for atmospheric bands, generalizing the square-root 
law. (21). 

We therefore approximate for small 0 



and 



1/sinh 0 = 1//3 - 0/6 
coth p = 1/0 + 0/3. 



By substituting into (26) and using the asymptotic 
expansion for J 0 , one obtains 



A = 1 - r - erf [(^/2)*] 

= erf [Vfl-wSa/d], 



with the usual definition 



erf 



( ) 2 f * 

V7T J 0 



e~ x2 dx. 



(27) 



(28) 



For small values of x, we have 

erf (x) = 2^/V^", 

and (27) reduces to the square-root absorption law 
(21). In practice, it is found that a square-root law 
holds only for a rather limited range of optical thick- 
ness ; then the curve for A bends around to reach unity 
asymptotically. It may be noted that in the derivation 



of (27) the region of a linear absorption law, which 
must always exist for sufficiently small u y has been 
obliterated. 

The bands of C0 2 consist of fairly regularly spaced 
lines, and for them the absorption law (27) is not too 
poor an approximation. On the other hand, the bands 
of water vapor and ozone have a highly irregular fine 
structure and cannot be well described by our approxi- 
mation. This model of the periodic band was intro- 
duced by Elsasser (1938, 1942) and is sometimes re- 
ferred to as the Elsasser band. 

4. Statistical band model 

We now treat another model for band absorption, 
based on statistical considerations, due to Goody 
(1952, 1954). Let Av be a spectral interval in which 
there are n lines of mean distance d\ 



Av = nd. 



(29) 



Let P(Si) be the probability that the ith line have an 
intensity Si and let P be normalized ; 



P(S)dS = 1. 



We assume that any line has equal probability of 
being anywhere in the interval Av. The mean absorp- 
tion clearly does not depend on v provided we are far 
enough away from the edges of the interval. We shall 
therefore determine it for v = 0, the center of the 
interval. We have from (14), if we let the center of the 
ith line be at vo — v iy 

ki(Si, Vi) = k{v = v i} VQ = 0). 

The mean transmissivity is found by averaging over 
all positions and all intensities of the lines ; thus, 

t = r dw - - r dv n 

(Av) n J to J A, 

r P{Sde-'"dSi- • • f P(S n )<r**«dS n . 
But all the integrals are alike; hence, 

= |\ - ^ J dv J P(5)(l - e-*«)dsj . 

By using (29), it is at once seen that the last expres- 
sion approaches an exponential when n becomes large ; 
we write 



r = exp 



[~^/ p(5)[ ^ iA ^4 (3o) 
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with the abbreviation 



DliA,]- f (1- 



(31) 



where A x is nothing but the absorptivity of a single 
line taken over the interval Av. The expression (30) is 
clearly valid for any line shape, k(S y v), whatever. 

We shall next consider two special cases. First, let 
all the lines have equal strength. The integral in (30) 
can simply be omitted and we have, by using 
again (29), 



(32) 



This reduces to the algebraic sum of the individual 
line absorptions, 1 — r = nA lt when the exponent 
is small, the lines absorbing weakly. Otherwise, the 
absorptivity is always less than nAi. This is due to 
the fact that overlapping of lines diminishes the ab- 
sorption and that in an irregular array there remain 
occasional wider absorption gaps which are filled up 
but slowly as the thickness of absorber increases. 

If we use the square-root absorption law (21), 
formula (32) becomes 



r = exp (— 2^uSa/d). 



(33) 



Next, let the lines be of different strength and 
consider a simple Poisson distribution for the prob- 
ability of their strength — namely, 



By letting k 



P(5) = — e -siso t 
So 

KS, we find from (30) 
KuSo 



(34) 



exp 



L d J 1 + KuSo J 



+ KuSo 

We now introduce the Lorentz shape (14)- 



K = 



t(v* + a 2 ) 



This vanishes fast enough for large | v j that we can 
extend the integration in (34) from — 00 to + 00 f 
giving 

' = ex P - - • (35) 

L dy! a 2 + uSoa/T J 



This is Goody's formula. It should be noted that it can 
be expressed as a function of two parameters only, 
apart from u, namely So/d and a/d, in complete 
analogy to the periodic band case (26). Usually, in 
atmospheric bands a « d, the lines are narrow com- 
pared to their mutual distance. If a is small, we may 
in a first approximation neglect a 2 under the square 
root, and (35) becomes 




Fig. 4. Comparing transmission for periodic band, equation (27) 
with random band, equation (36) both for limit of widely spaced 
lines. (0 is the error integral (28).) 



which differs from the corresponding formula (33) 
for lines of constant intensity only by a slight numeri- 
cal factor in the exponent. 

In fig. 4, the two functions (27) and (36) are com- 
pared, with a logarithmic scale for the argument, 
which is the same for both cases provided we identify 
5 in the former with So in the latter. Actual band 
spectra are neither as regular as is implied in the one 
formula nor as irregular as is assumed in the other; 
hence, one might expect most real transmissivity 
functions to lie between these two cases. 

The basic formulas (30) and (31) of the statistical 
transmission theory are very general, and it is possible 
to deal with sets of lines having a shape different from 
the Lorentz shape. The Doppler shape becomes im- 
portant in the upper layers of the stratosphere. This 
will be dealt with in section III. Calculations of this 
more general type have been carried out by Kaplan 
(1953) who has also used this method to treat the 
statistical superposition of several period bands of the 
type described in section II ; see, also, Plass (1958). 

5. Pressure. Semi-empirical method 

We shall here merely state the facts regarding the 
pressure effect on absorption; some of the physical 
background will be discussed in section III. The 
finite width of a Lorentz line is produced by the col- 
lisions of the radiating molecule with other molecules ; 
from kinetic theory, the number of collisions is pro- 
portional to p/ifT; we write, therefore, 



P \To 

Po * r 



(37) 



r = exp (— -yJiruSoa/d), 



(36) 



where the width ao refers to standard conditions, po, 
To. Since the relative concentration of the optically 
active gases in the atmosphere is usually quite small, 
p essentially designates pressure of the ambient air. 
The pressure effect expressed by (37) is very important 
in atmospheric transfer calculations ; the temperature 



8 



METEOROLOGICAL MONOGRAPHS 



Vol. 4, No. 23 



effect in (37) is much smaller and may be ignored 
when no great accuracy is required. 

Let us now look at a relatively simple transmissivity 
formula such as (35). Although we have succeeded in 
smoothing over the rapid variation of absorption 
caused by the numerous lines, there remains a de- 
pendence on some parameters, and the application of 
such a formula to a great variety of atmospheric 
conditions will be too cumbersome to be practical. 
On the other hand, much simplification occurs when 
the lines are very narrow. Then (36) applies for the 
statistical band and (27) for the periodic band. Both 
depend only on one combination of parameters. 
In the atmosphere, the spectral lines become the 
narrower the higher we go, as seen from (37) ; since 
the lines are already fairly narrow at one atmosphere 
of pressure, they will be even more so higher up in 
the atmosphere. This suggests that, in a fair approxi- 
mation, transmissivity functions may be expressed as 
generalizations of (36), say 

t — r(u*L) f 

where L is designated as generalized absorption coeffi- 
cient and is a slowly varying function of v over the 
band. In order to take account of the pressure de- 
pendence as given by (37), we note that only the 
product ua appears in (36) and (27). We therefore 
introduce a "reduced" optical thickness, 

M * = M fv?' (38) 

and write 

r = T (u*-L). (39) 

It is possible to check the accuracy of (39) by means 
of laboratory measurements. Numerous wide-slit ob- 
servations on the relevant far infrared bands are now 
available. These are made with variable optical path 
lengths, u t and with variable air pressures. By fixing on 
a given interval of the spectrum, one can then empiri- 
cally obtain the function (39) and can check how far 
the pressure reduction (38) represents the actual data. 
By doing the same thing for a number of spectral 
intervals, one can furthermore see how much the 
different r's for various intervals differ from each 
other. It is then possible to justify the use of a single 
r for the spectrum of one substance, as it appears that 
the different r's for different intervals are rather close 
to each other. Once this r is adopted, one can de- 
termine the variation of L with frequency over the 
band. This is simply done by sub-dividing the band 
into a number of intervals (much smaller than the 
total band width but still wide enough so that the 
individual lines are completely smoothed out). One 
then determines L for each interval so that the result- 
ing functions (39) fall upon each other as closely as 



possible; one next draws a smooth curve through 
these L-values as function of v. This is the desired 
curve of the generalized absorption coefficient. 

By means of the laboratory observations, one can 
check directly the accuracy of this procedure. Details 
for particular gases will be given in section IV. We 
should note here already that the appropriate scale is 
logarithmic throughout. Thus, we shall give (39) 
the form 

r « r(logw* + logL), (40) 

and we shall use decadic logarithms in all cases. The 
absorption is found to be widely distributed over 
several cycles of the logarithmic scale, and it would be 
impossible to do justice to the observations with a 
linear scale. As a result of this, no great accuracy is 
required in the determination of the optical length, u. 
This is advantageous in applications to the actual 
atmosphere, since it means that no great accuracy is 
required in determining the distribution of the opti- 
cally active gases as function of height. 

6. Stratified medium 

So far we have assumed straight beams of radiation 
passing through an absorbing medium. This corre- 
sponds to laboratory experiments with absorption 
tubes but does not of course represent the geometry of 
the atmosphere. We assume the atmosphere to be 
horizontally stratified and of infinite extent in the 
horizontal direction. Let the axis of the polar coordi- 
nates <p be in the vertical. It is convenient to inte- 
grate the radiative currents (beams) over a hemi- 
sphere. This leaves us with two integrated currents, 
one whose vertical component is downward and one 
whose vertical component is upward. These currents 
integrated over the hemisphere will be designated as 
fluxes ; thus, there are an upward flux and a downward 
flux, both obtained by integrating straight-line beams. 

Given a level in the atmosphere, it will be crossed 
by a family of beams with upward components of in- 
tensity /(#), assuming independence oi<p. Simple 
geometry then shows that the net upward flux is 

/■ir/2 

/ = 2tt I /(*) cos 0 sin (41) 
Jo 

and similarly for the downward flux. For black-body 
radiation, / is independent of and we find 

f b b B, = irh, (42) 

which is identical with (3). 

Next consider a layer of optical thickness u upon 
which isotropic flux io, (dIo/d& = 0), is incident on 
one side. If there was only a vertical beam, the inten- 
sity emerging on the other side would be / = hr{u). 
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By integrating now over the hemisphere, we see with 
the help of fig. 5 that 

t{u sec &) cos $ sin &d&. 

o 

We write this, by using (42), as 

/ = Mriu), (43) 

where tf will be designated as the flux transmissivity 
and is defined by 

n*ft 

t f (u) = 2 | r(u sect?) cos # sin &d&. (44) 

We can rewrite this as 

r r (u) = 2u* j r(x)dx/x\ (45) 

which lends itself rather readily to numerical com- 
putations. Since the transfer equations used here are 
linear, the integration over the angles may be carried 
out at any stage of the integration, and hence flux 
transmission in the atmosphere can be taken into 
account by a simple formal device — namely, by 
treating the atmosphere as if there were only two 
beams of radiation, one going vertically upwards, the 
other vertically downwards. It suffices to replace in 
the beam formalism / by / and, at the same time, r 
by tf. With these substitutions, the transfer formulas 
for the two vertical beams carry over unchanged into 
those for the two vertical fluxes, as may readily be 
proved explicitly. 

The shape of tf is in all practical cases extremely 
similar to that of the beam transmissivity, r, except 
that there is a slight shift on the log u scale. For most 
practical purposes, it suffices to set 

tf(u) = r(1.60w) 

or, equivalently, 

TF(log u) = t (log u + 0.20). (46) 

In the work described below, we have usually carried 
out the integration (45) for the empirically obtained 
beam transmissivities of the atmospheric gases, but 
the advantage gained thereby over (46) is small. 

7. Uniform isothermal layer 

The following results are due to Pederson (1942) 
and, independently, to Strong and Plass (1950). 
Although they are not used in the present work, we 
quote them here for the sake of completeness. They 
refer to the case where the integration over the height 
in the atmosphere can be carried out analytically. This 
can be done in a model, representing the stratosphere, 




Figure 5. 



where an isothermal layer of uniform chemical com- 
position extends from a lower boundary (indicated by 
the subscript 0) upward to infinity. For an isothermal 
layer, (37) reduces to 

a = ctop/po. (47) 

Within the layer, the density decreases upward accord- 
ing to the barometric law. Since there ir> uniform 
mixing, we have, for the optical thickness in the 
vertical, w, 

du = cpdz 

where c is the mixing ratio of the optically active gas 
and p the density of the air. Hence, by the barometric 
law, 

gdu — cdp = cpoda/ao, (48) 

where g is the acceleration of gravity. We compute the 
absorptivity of a single line, by using (14). We slightly 
generalize by admitting oblique beams making an 
angle # with the vertical. The absorptivity of the layer 
is found after some calculations to be 

2a f I f x 2 + (ao/a) 2 V 1 

am -tJ. I'-l-JTi-Jr (49) 

where 

7 = Scpo sec$/2-irgaQ. (50) 

If the term (ao/a) in the numerator is neglected (corre- 
sponding to wing absorption), the integral can be 
simply evaluated by means of gamma functions. Plass 
and Fivel (1955) have, moreover, shown that (49) can 
be expressed in terms of Legendre functions — namely, 

27T7a 0 y 

A U*o, a) = 

X [p t (cosh In — ^ P y _, (cosh In — ^ J. 

As Goody (1954) has pointed out, such expressions 
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are of limited applicability because they do not take introduced the periodic Elsasser band into this model 
into account overlapping; and, indeed, A goes to of the stratosphere and succeeded in analytically 
infinity for infinite y. The same authors have, however, evaluating the resulting integrals. 



III. PHYSICS OF THE ATMOSPHERIC BANDS 



1. Spectroscopic characteristics 

This section is intended to justify in more detail the 
techniques and approximations used from the view- 
point of molecular physics. These arguments mainly 
corroborate the validity of the formalism already pre- 
sented, and most of it may therefore be passed over 
rapidly by those interested only in the meteorological 
aspects of the problem. However, section 1 1 1.4 deals 
with the modification of the band-absorption contour 
by the temperature which will be used extensively 
later on. 

The optically active gases of the atmosphere, C0 2 , 
H 2 0, and O3 are all triatomic molecules. The band- 
spectroscopic properties of such molecules are des- 
cribed in the well-known work by Herzberg (1945) to 
which we may refer the reader for all particulars. We 
shall here give just such details as are desirable for the 
consistency of our own account. 

Carbon dioxide. Spectroscopic evidence indicates 
unambiguously that the three atoms of C0 2 form a 
symmetrical straight-line array, the carbon atom in 
the middle being flanked by oxygen atoms on either 
side. Such a configuration cannot have a static electric 
dipole moment. An electric dipole moment will, how- 
ever, be induced by vibrational motion of an asym- 
metrical type. Fig. 6 shows the three normal modes of 
vibration of such a configuration. From what has been 
said, it follows that the symmetrical motion vi should 
not give rise to an electric dipole moment and should 
not, therefore, be optically active. It has been iden- 
tified in the Raman spectrum near 7.5 /x. In the vibra- 
tion V2, the electric moment is perpendicular to the 
axis of the molecule. The 15-ju band represents this 
particular vibration, the band being a "fundamental," 
that is a transition from the ground state to the first 
excited vibrational state. Another fundamental corre- 
sponds to the transition of the first vibrational state, 
v 3 . This is the 4.3-/z band which appears at the short- 
wave edge of the black-body curve for atmospheric 
temperatures. 

A vibrational level diagram of the C0 2 molecule may 
be found in Herzberg's book. The only transitions of 
interest to us here besides the one already mentioned 
are some rather faint bands which appear near 10 n 



and which, in spite of their weakness, might have a 
not entirely negligible effect upon heat transfer in the 
atmosphere. Band-spectroscopically, they are vibra- 
tional harmonics. We shall revert to the details in 
section IV. 

Owing to its straight-line arrangement, the C0 2 
molecule has spectroscopics closely akin to those of a 
diatomic molecule. The 15-/* band in particular has all 
three branches of a typical band of such a molecule, 
the P, Q, and R branch, corresponding to a jump of 
the rotational quantum number by +1, 0, and — 1. 
In a very crude approximation, the lines of the P and 
R branches are equidistant ; thus, if it was not for the 
superposition of these two branches, the periodic 
Elsasser model would be an excellent approximation 
to the actual transmissivity. As it is, the observed 
transmissivities are somewhat changed toward the 
random Goody model. The Q branch is quite closely 
clustered near the center of the 15-/* band and for that 
reason does not contribute importantly to radiative 
transfer. It will not be taken into account separately 
in the semi-empirical developments of section IV. 

Water vapor. The water molecule forms an isosceles 
triangle which is obtuse and has an apical angle of 
104.5 deg. Fig. 7 shows the three normal modes of 
vibration for such a structure. The 6.3 band has been 
identified as the fundamental v 2 of fig. 7. The two 
fundamentals, vi, and *> 3 , are found close together, in 
a band near 2.7 ju, well on the short-wave side of the 



V, ♦ v 2 ♦ v 3 

Fig. 6. Normal modes of a linear molecule. 




Fig. 7. Normal modes of molecule forming an isosceles triangle. 
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Fig. 8. A section of the rotational water- vapor band after Randall, Dennison, Ginzburg, and Weber (1937) taken with large 
spectral resolution. Numbers are cm" 1 ; triangles above spectrum give lines computed from quantum mechanics by position and in- 
tensity ( = area of triangle). 



region which interests us here. This band and the 
bands representing higher harmonics are of meteorolo- 
gical interest in that they give rise to absorption of 
sunlight by the atmosphere. 

The band on the long-wave side of 20 n represents 
the purely rotational spectrum of water vapor. The 
water molecule forms, with respect to rotation, an 
asymmetrical top, and the line structure of the spec- 
trum does not have the simplicity of a symmetrical 
rotator which we found in the CO2 molecule. There is 
no clearcut regularity at first sight ; a typical stretch 
of the rotational band as it appears with very high 
spectroscopic resolution is shown in fig. 8 (Randall 
et al., 1937). Above each section of the spectrum, there 
is a line studded with triangles which represent posi- 
tions and intensities ( = areas of triangles) calculated 
from quantum mechanics. Such a spectrum may be 
approximated by the Goody transmission formula, 
and the laboratory observations indicate that, with 
suitably chosen parameters, formula (35) can be made 
to give an excellent fit to observed transmissivities 



taken with a sufficiently wide slit. The fine structure 
of the 6.3-fx band is essentially similar to that of the 
pure rotational band. 

There is no clear evidence that water vapor has 
absorption bands in the region between the two bands 
mentioned — that is, between about 8 n and 13 /*. We 
shall return to this point in section IV. 

Ozone. The analysis of the far infrared spectrum of 
ozone remained unsatisfactory for many years, except 
that it seemed certain that the molecule forms an 
isosceles triangle. In 1948, Wilson and Badger rein- 
vestigated the spectrum and discovered a faint band 
centered at 9.0 hitherto overlooked. If this is inter- 
preted as the fundamental v x of fig. 7, then the well- 
known absorption peak at 9.6 n can without difficulty 
be interpreted as y 3 , and the 14-ju band, which is some- 
what weaker than the latter, can be interpreted as v 2 . 
These assignments are corroborated by the more recent 
investigations of Hughes (1952-53) who finds an 
apical angle of 117 deg and by the analysis of Kaplan, 
Migeotte, and Neven (1956). The infrared bands of 
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ozone are of interest for the radiative cooling of the 
stratosphere where their thermal effects may be quite 
important. In our semi-empirical approach of the next 
section, we shall deal in detail with absorption in the 
9- to 11-ju region. We shall merely give absorption 
data for the 14-/J band without entering into the 
problem of overlapping with the C0 2 band at the same 
location. Likely, no appreciable error is committed in 
most cases if the 14-/* band of ozone is omitted alto- 
gether in radiative-transfer calculations. 

2. The shape of spectral lines 

The Lorentz shape (14) of the far infrared lines is 
fundamental for the whole theory of atmospheric 
transmission. We shall therefore give an approximate 
treatment, showing how this distribution of intensity 
in a spectral line arises (see, also, Born, 1933, p. 435). 
An isolated molecule would emit an almost purely 
harmonic wave, say Ae 2 * ivt y and the only deviation 
from purely harmonic behavior would be produced by 
the damping that results from the emitted wave itself ; 
in the far infrared, the spectroscopic effect of this 
damping is exceedingly small. During the time a 
radiating molecule is in collision with another molecule, 
the radiating wave train is altered ; owing to the inter- 
molecular forces, the frequency of the emitting mol- 
ecule is temporarily shifted by an appreciable amount. 
Now the fraction of time a molecule spends in collision 
with other molecules is very small (about 10~ 3 at STP). 
Lorentz therefore assumes that the principal effect of 
the collision is to destroy the phase coherence of the 
emitted wave train ; one may assume the collision to 
be instantaneous, but such that after the collision the 
molecule starts emitting with another phase, the new 
phases being randomly distributed. This is equivalent 
to saying that the spectrum of an assembly of mol- 
ecules consists of a set of wave trains of variable 
lengths with arbitrary initial phases. From general 
statistical principles, the time between collisions is 
distributed according to an exponential (Poisson) law, 
the probability that a collision occurs between / and 
being e~ tltQ where U is the mean time between 
collisions. We must also average over all the initial 
phases of the wave trains, and we can express this by 
using not the spectrum itself but the power spectrum 
(the absolute value squared of the Fourier transform). 
Hence, the desired spectrum is 



If we set l//o = 2wa, this becomes equivalent to (14) 
apart from constant factors and apart from the de- 
pendence on 1/p in front of the right-hand side of (51). 
However, since v varies but very little over the line 
width, this factor may be set equal to a constant in an 
excellent approximation. Now, since 2-kol is the number 
of collisions per molecule per unit time, the dependence 
of a on pressure and temperature is, from the kinetic 
theory of gases, just given by (37). 

An elaborate quantum-mechanical theory of Lorentz 
line broadening has been given by Van Vleck and 
Weisskopf (1945). Their calculations apply more 
closely to lines in the visible, rather than in the far, 
infrared. An extension of the theory to the latter 
region was given by Anderson (1949) who showed that, 
in addition to broadening of the line, pressure produces 
a slight shift of the line center. Plass and Warner 
(1952) have investigated the effect which this has on 
atmospheric transfer and find it as a rule negligible. 
If the time which the radiating molecule spends in 
collision is not completely negligible, there appear 
deviations from the Lorentz shape at large distances 
from the line center. References may be found in the 
paper by Plass and Warner where it is also shown that 
these latter effects may be neglected in atmospheric 
transfer. Anderson found that the line width varies 
to some extent as a function of the rotational quantum 
number in each sub-band and also from sub-band to 
sub-band. In our preceding theoretical developments, 
it has been assumed that all lines have the same width. 
It seems safe to say that such added analytical com- 
plexity as would be introduced by accounting for a 
width varying from line to line would not be repaid by 
significantly novel physical results. We confine our- 
selves to interpreting the formal line width, a, of our 
transmission theory as the average width of the lines. 

A second cause of line broadening next to collisions 
is the Doppler effect owing to the motions of the 
radiating molecules. Let us assume for the moment a 
highly rarefied gas where there is no collision broaden- 
ing. All molecules in a given quantum state radiate 
at wave number ? 0 , but, as the theory of relativity 
shows, this appears shifted as seen by a stationary 
observer to the wave number 
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v = fo(l db v/c), 



(52) 



where v is the velocity of the molecule relative to the 
observer. Now the fraction of molecules between v 
and v + dv is, from kinetic theory, 



n(v)dv — 



2kT 



e ~ mv lj2kT ^ 



(53) 



(51) 



where the symbols have their usual meaning. To obtain 
the Doppler distribution, we merely replace in (53) 
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the variable v by v, by using (52); if, moreover, we 
normalize to an integrated line intensity 5, we readily 
find for the distribution of the absorption coefficient 

™ i^Yl <s4) 

where 



c m 

is a measure of the Doppler width of the line, the half- 
width at half maximum being 5(/w2)*. 

In the upper atmosphere, we find a combination of 
collision broadening and Doppler broadening, whereas, 
in the lower atmosphere, Doppler broadening is 
usually much less than collision broadening and may 
safely be neglected in the computations. When both 
broadening effects are present simultaneously, the 
resulting line shape is a combination (mathematically 
speaking a convolution) of the two (see Mitchell and 
Zemansky, 1934). Let K be a shape function nor- 
malized to unity, so that 

kh = SKl> kD — SKd 

where kh and Ub are defined by (14) and (54) respec- 
tively ; then the resulting line shape of a line of total 
intensity S will be 

k(v - fo) = S j dv'K L (v - v')K D (v' - *o). 

This function has been exhaustively studied by Plass 
and Fivel (1953). Their findings may be summarized 
by saying that, on the whole, the effect of Doppler 
broadening upon a Lorentz-broadened line is small. 
We may distinguish here between "strong" and 
"weak" lines. This terminology refers to the Strong- 
Plass formalism discussed in section 1 1.7 which is 
more realistic for the atmosphere than a uniform 
absorbing layer. A line is strong when y » 1 and weak 
when y « 1, where the parameter y is defined in (50). 
The intuitive physical meaning is as follows: in the 
Strong-Plass atmosphere, the line gets narrower as 
one goes upwards, owing to the decrease in half -width 
proportional to the decrease in p. This produces an 
intensification of the absorption in the line center 
which makes up for the decrease of optical thickness 
in the upper layers. For a weak line, the absorption 
remains always in the region of the linear law (20). 
For a strong line, on the other hand, the center is 
always completely opaque, and absorption occurs in 
the wings so that the square-root absorption law (21) 
holds. For a weak line, it may be shown that the total 
absorption is independent of the mechanism of 
broadening, a plausible enough result. In atmospheric 
transfer, the strong lines (and, in particular, the very 



strong lines) are much more important. Now, as we 
go up in the stratosphere, the Lorentz width decreases, 
whereas the Doppler width, which depends only on 
VT, is fairly constant. For any line, there exists a 
level at which the Lorentz width and Doppler width 
become equal. Now, if the line is strong, transfer takes 
place in its wing, whereas no transfer occurs in its 
central parts. The intensity of a Lorentz line decreases 
only slowly, as (p — ? 0 )~ 2 , with distance from the line 
center, whereas the Doppler line decreases exponen- 
tially. Thus, the effect of Doppler broadening in the 
wings should be small for strong lines. Since most of 
the lines contributing appreciably to atmospheric 
transfer are strong, this means that the approximation 
in which Doppler broadening is neglected should still 
be very good, well above the height at which the two 
widths are equal. This is quantitatively confirmed by 
the calculations of Plass and Fivel. 

So far, we have only dealt with pure theory. The 
measurement of the detailed intensity distribution of 
a line in the far infrared is exceedingly difficult, at the 
limits of the possible, because the only spectral sources 
available are thermal ones, and their intensity is not 
far enough above noise to permit unlimited resolution. 
Numerous experimental investigations of line broaden- 
ing have been carried out, however, in the visible and 
in the near infrared, and the Lorentz distribution (14) 
has been well verified under the physical conditions 
where it applies. It has also been verified in the micro- 
wave region where good resolution is readily available. 
In the far infrared region, it has been verified so far 
as feasible ; the line width in the far infrared atmos- 
pheric bands is of order 0.05 to 0.15 cm -1 at p = 1 atm. 
The Doppler width, on the other hand, is much 
smaller; it is several times 10~ 3 cm" 1 at atmospheric 
temperatures. Goody (1954, p. 156) has gathered 
numerical data on line width (see also Kaplan and 
Eggers, 1956). Fortunately, the precise values of the 
width of individual lines are not significant for our 
semi-empirical approach. The height in the atmos- 
phere up to which Doppler broadening may be 
neglected as compared to Lorentz broadening is of 
importance. Goody gives the following heights at 
which a = 5, the Lorentz width equals the Doppler 
width, for various atmospheric bands: 23 km for the 
6.3-/* water-vapor band, 35 km for the rotational 
water-vapor band, and 34 km for the C0 2 band at 15 ju- 
Furthermore, the work of Plass and Fivel (1953) indi- 
cates that in the case of the stronger lines the Doppler 
width may be neglected even somewhat beyond the 
point where it becomes equal to the Lorentz width. 
In the case of very strong lines (7 > 100) which carry 
an appreciable part of atmospheric transfer, Plass and 
Fivel calculate that Doppler broadening should be 
negligible up to a height of 50 km. In view of the above 
figures, it seems safe to say that taking into account 
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Lorentz broadening only and ignoring Doppler 
broadening should yield a reasonable approximation 
up to heights approaching 40 km. As a matter of fact, 
the determination of radiative transfer in the upper 
stratosphere does not depend so much on the spectro- 
scopic complexities just mentioned ; the main problem 
is that the quantities of optically active gas present 
at these altitudes are so minute that it becomes nearly 
impossible to reproduce these conditions in the 
laboratory. In other words, the integrated curves of 
total band absorption which will be given later are 
likely to be rather inaccurate at their lower end. This 
is likely to be a much larger source of error than any 
other. 

3. Local equilibrium and blackness 

We shall next discuss the question as to how far 
Kirchhoffs law (1) is quantitatively valid for far 
infrared atmospheric radiation. The transfer equations 
used here are based entirely upon this law. In most 
radiative-transfer problems involving visible light, as 
encountered in astrophysics as well as in certain 
technical applications, there is an appreciable amount 
of scattering of the radiation. This leads to a transfer 
formalism much more complicated than the one used 
here (Chandrasekhar, 1950). In the far infrared, on 
the other hand, scattering is unimportant except 
perhaps for solid particles. We shall not enter here into 
a discussion of the detailed mechanisms of scattering 
for which quantum mechanical theories are available 
(Heitler, 1954). The scattering which occurs in the 
atmosphere cannot be attributed to individual mol- 
ecules but is due to the statistical density fluctuations 
of the air molecules (Rayleigh Scattering ; see Born's 
Optik, p. 371 ff.). In any event, molecular scattering, 
whether individual or collective, is inversely propor- 
tional to the fourth power of the wavelength. This, as 
is well known, is the reason that the sky is blue, 
scattering being important at the short-wave end of 
the visible spectrum and becoming rather weak at its 
long- wave end. As we go to the infrared, it becomes 
still weaker, and, in the far infrared, molecular scatter- 
ing is totally negligible, as simple calculations by the 
Rayleigh formula show. 

Next in importance to the absence of scattering is 
the existence of local thermodynamical equilibrium. 
By this is meant the fact that the thermodynamical 
state of the atmospheric gas as a whole (air plus 
optically active substances) may be characterized by 
a single temperature, T, so that the mean population 
of any state of energy E is given by the Boltzmann 
formula, n(E) = P _1 exp( — E/kT) f where P is the 
partition function. In order that this be so, it is 
necessary that when the population of any particular 
quantum state or group of quantum states becomes 



enriched or depleted, this change is rapidly distributed 
over all the quantum states of the whole gas so that 
no accumulation in any one state occurs. In the con- 
crete circumstances of infrared radiation transfer, this 
means that the molecular levels which become excited 
by the absorption of the infrared flux must redistribute 
their excitational energy by collisions before they 
re-radiate it. This means that the relaxation time of 
rotation-vibrational energy must be short compared 
to the lifetime of the excited molecular states. A close 
analysis of the available data shows this to be the 
case. The quantitative information has been discussed 
extensively by Goody (1954, p. 146), and we shall 
therefore not repeat the detailed arguments here. 
Goody estimates that local thermodynamical equili- 
brium prevails in the atmosphere up to a level of 50 
km. It may therefore safely be assumed that within 
the range of validity of the present tables any devia- 
tions from local thermodynamical equilibrium will be 
far smaller than many other errors involved. 

In dealing with the solution of radiative- transfer 
problems, one must, as a rule, fulfill boundary condi- 
tions. The most important boundary is of course the 
ground, but the top and bottom of a cloud layer are 
also significant. Experience shows that in the fai 
infrared the optical properties of most solids or liquids 
become extremely simple; most solids are essentially 
black, excepting only the metals and certain crystals. 
By a black body, we mean again one which absorbs 
completely all radiation failing upon it, independent 
of frequency. Conventional ground is black to within 
the limits of error of most measurements. Dorsey 
(1940) gives values of reflectivity of liquid water as 
function of wave length at normal incidence : approx- 
imately 2 per cent at 5 /x, 1 per cent at 10 jjl, 5 per cent 
at 15 ja, and 10 per cent at 60 ix. So far, the calculations 
of atmospheric transfer have always assumed that the 
lower boundary of the atmosphere is black. Also, the 
errors committed in assuming the same for the sea 
surface appear small enough so that considering the 
sea as a black body is usually a good approximation. 
To justify this further, note that opaque regions of the 
atmosphere are those for which the free path of radia- 
tion between emission and absorption is small. Radia- 
tive transfer is then also small, and some reflectivity 
of the water surface is unimportant. On the other 
hand, in the regions of the spectrum where the moist 
atmosphere is relatively transparent, the radiation 
travels a long way, and radiative cooling becomes 
important. The measurements quoted by Dorsey show 
that the reflectivity of liquid water is large in the 
regions where the water-vapor phase has absorption 
bands but small in the region centered about 10 ju 
where vapor is transparent. This argument would 
seem to indicate that the effective reflectivity of the 
sea surface in radiative-transfer problems is much 
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smaller than would be indicated by a mere arithmetical 
average of the above figures for reflectivity. 

Observation shows that a cloud layer (other than 
cirrus) which is more than a few meters thick is com- 
pletely black in the far infrared. Clouds, therefore, are 
considered in the computations as black boundaries 
of the radiative fluxes. Since the top of a cloud as a 
rule cools rapidly by radiation into space, the lapse 
rate in a cloud is usually adiabatic, and heat transport 
takes place by convection. 

4. Temperature dependence of band contour 

In section 1 1.5, we introduced the generalized 
absorption coefficient L(v) such that the transmis- 
sivity as function of wave length and optical path is 
given by t{u*L) where w* is the pressure-reduced 
optical path. If it be assumed that L is independent of 
temperature, then the temperature appears in the 
transfer equations only via the black-body flux B v 
which of course is strongly temperature dependent. 
There is, however, a secondary dependence on tempera- 
ture, through L(v) itself. We shall explain how this 
temperature effect arises and how it is taken into 
account in our computations. 

The change of absorption with temperature is due 
to the change of the population of the energy levels of 
the molecules with temperature. The population of 
the ith level, of energy E i} is given by the Boltzmann 
formula. 

»(£#) = P^ e ~ E ^ kT (56) 

where P is the partition function, the sum of the 
Boltzmann factors of all states. Now the excited 
rotational states of a molecule at atmospheric tem- 
peratures are very numerous, as evident from the 
large number of lines in a band. Short of carrying out 
extremely elaborate quantum mechanical calculations 
which would take into account all these levels, only 
a rather schematical model of molecular behavior can 
provide us with the necessary correction terms in a 
simple enough form. 

The model adopted here is as follows. Assume for a 
moment that we replace the band spectrum by the 
Fourier spectrum of a symmetrical rotator in classical 
physics. This latter spectrum is obtained by distribu- 
ting an ensemble of classical symmetrical rotators over 
the frequency range involved according to a Boltzmann 
formula. The kinetic energy of such a rotator is 

E = (2tt^) 2 //2, 

where I is the amount on inertia. Thus, the Boltzmann 
distribution of such an ensemble is given by 

n(E) = i«exp(- cvVT) (57) 



where c is some constant. This shows that if we plot 
log n against v we get an inverted parabola. In this 
simple model, the generalized absorption coefficient 
would be 

L{v) = const. v z n{v). (58) 

The variation with *> 3 arises from radiation theory for 
fixed n ; it is, however, much less rapid than the expo- 
nential variation of n(v) ; hence, in the model log L(v) 
appears in a first approximation simply as an inverted 
parabola. Now the two sides of a rotation-vibration 
band are generated by adding and subtracting, 
respectively, the frequency of rotation to the frequency 
i>o of the vibrational jump. Hence, we might ask 
ourselves whether \ogL(v) of a molecular band may 
be schematized by means of a full inverted parabola. 
It turns out that this is somewhat too crude an approxi- 
mation, mostly owing to the asymmetry of real bands 
(produced by centrifugal stretching of the molecule). 
We have therefore adopted a scheme in which the two 
sides of a band are approximated by different parabolas 
(there is of course only one "side" to the rotational 
band of H 2 0). We thus let 

logL = - a(v - vo) 2 + b (59) 

where we have neglected the dependence of L on v z 
expressed in (58) as compared to the exponential 
dependence of n{v) on v. Here, vo is the band center. 
The parameters a and b were obtained for each side 
of a band from the measured data (section IV) by a 
graphical triai-and -error method. This sufficed since 
the fit of the parabolas to the actual contour of log L 
in any event is only moderately good. On the other 
hand, it should be remembered that the parabolas 
themselves are not used but only certain differences 
derived from them: the band contours observed at 
room temperature are corrected by means of the 
temperature dependence of these substitute parabolas. 

To find the temperature dependence, we notice first 
that by the definition of P we have from (56) for a 
simple rotator 

P = f°exp(- E/kT)dE = kT. 
Jo 

In view of this, we use (57) and (58) to define the 
variation of the generalized absorption coefficient on 
each side of a band, neglecting again the dependence 
on v z in (58). Thus, 

const. f c 1 
L{v - v 0 ) = ~y~ exp [ - - (v - »o) 2 J (60) 

where c is a constant. Now let L 0 refer to laboratory 
temperature r 0 (293K, say). Then (60) readily gives 
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the desired temperature-correction formula — namely, 

T Q - T To 
log L = log L Q - a — (v - vq) 2 + log — , (61) 

where a — c/T 0 . Moreover, if in (59) we let L = L 0 , 



the present a becomes identical with that of (59). 
Thus, on obtaining a and b as used in (59) by a trial- 
and-error process, the correction (61) is fully deter- 
mined. The method outlined was applied numerically 
to all the bands for which computations were made. 



IV. LABORATORY TRANSMISSIVITIES 



1. Formulary 

The quantities to be determined from laboratory 
measurements are the transmissivity, r, and the 
generalized absorption coefficient, L. By (39.) r is a 
function of the product, u*L f where u* is the reduced 
optical path given by (38). In the laboratory, the 
temperature is effectively constant because it is not 
feasible to lower the temperature in the rather bulky 
absorption apparatus. Hence, it suffices to write 

u* = up/po, (62) 

where p is the total pressure (almost entirely air 
pressure as a rule) and po is the standard pressure of 
76 cm Hg. In the remainder of this monograph, we shall 
for simplicity omit the asterisk and write ufor the reduced 
path throughout. 

We shall consistently use a logarithmic scale for u 
as well as for L and therefore write 

r = r(log^ + logL). (63) 

All logarithms in this work are to base 10. In the graphs 
of this section, where r is plotted as a function of its 
logarithmic argument, we simply use log u as the 
abscissa, it being understood of course that this 
designates log u + log L. 

The unit of u has been chosen differently for C0 2 
and 0 3 on the one hand and for water vapor on the 
other. For carbon dioxide and ozone, u is measured 
by the length in cm of a column of the pure gas at 
standard temperature and pressure. For water vapor, 
u is measured in cm of precipitable water — that is, 
grams per sq cm of the column. The latter unit is of 
course chosen so as to conform with common meteoro- 
logical practice. 

To determine r for any given region in the spectrum, 
one chooses an interval wide enough to comprise a 
large number of individual spectral lines (e.g., by 
taking the slit of the spectrograph wide enough) and 
plots the observed mean transmission of the interval 
against log u for a variety of values of u y where u now 
is the path length reduced for pressure according to 
(62). If this is done for a number of intervals in the 
spectrum, it is found for the gases considered here that 
the curve has closely similar shapes for different in- 



tervals of one and the same gas. One then adopts a 
standard transmissivity curve for any one gas. This 
curve undergoes a shift along the log u scale when one 
passes from one interval in the spectrum to an adjacent 
interval. If this shift is plotted as a function of v, one 
obtains a series of points which thereafter are replaced 
by a smooth curve approximating them. This is the 
curve of log L where L is the generalized absorption 
coefficient. 

The functions r(u) and logL(i>) provide the mate- 
rial wherewith to construct a radiation chart as des- 
cribed in the next section. Different substances have 
different transmissivity curves ; most empirical trans- 
missivities can be fitted by a Goody curve of type (35) 
on suitable choice of the two available parameters; 
this is particularly true of H 2 0 and 0 3 . Again, C0 2 
has a somewhat simpler band structure, and its trans- 
missivity curve does not differ too widely from the 
simple error integral (27). In the case of water vapor, 
we have used the same transmissivity curve for both 
bands, the 6.3-/* band as well as the rotational band, 
but the region between the two bands, the "window/' 
requires a more involved treatment. 

With the exception of the values for the water- vapor 
window which are in part the result of computations, 
all our curves are based on laboratory measurements 
only, as made with equipment allowing long optical 
paths and the addition of air to study pressure effects. 
For carbon dioxide and water vapor, the measurements 
have been made at two laboratories, and, in fact, inde- 
pendent sets of data are available in the case of each 
band so as to provide a mutual check. These are the 
Laboratory of Astrophysics and Physical Meteorology 
at Johns Hopkins University under the direction of 
Professor John Strong and the Physics Department of 
Ohio State University with its old tradition of infrared 
spectroscopy. The measurements for ozone come from 
different laboratories ; there are also two independent 
sets. 

2. Carbon dioxide 

The experimental work on the 15-m carbon-dioxide 
band used here has been available only in the form of 
technical reports. A report by Cloud (1952) was used 
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Fig. 9. The 15-/* carbon-dioxide band taken with a wide slit. Type of data used for our computations. 



by Elsasser and King (1953) to derive a transmissivity 
curve and a curve of the generalized absorption coeffi- 
cient which are presented below. The first computation 
of the C0 2 radiation chart was based on these values. 
It was ultimately modified by using further, independ- 
ent, experimental data contained in a report by 
Howard, Burch, and Williams (1955) and by 
Yamamoto and Sasamori (1958). 

The work of Cloud was done with the 100-ft trans- 
mission tube at Johns Hopkins University. Mirrors 
at the ends of the tube provide a six- fold traversal, the 
total path length being 190.5 m. Nitrogen was admit- 
ted as a line-broadening agent, which may be shown 
to have for all practical purposes the same broadening 
effect as atmospheric air. In some cases, helium was 
used instead; it gives very nearly the same effect as 
nitrogen. Both the partial pressure of C0 2 and that of 
were varied between extremely wide limits, cover- 
ing almost all the conditions to be met with in the 
atmosphere. The 15-/* band was scanned with a very 
wide slit which obliterated all finer detail of the spec- 
troscopic structure. A series of traces obtained by 
Cloud is shown in fig. 9. The C0 2 pressure in this 
particular series was 0.2 cm of Hg. The partial pres- 
sures of nitrogen and helium were as indicated for 
each diagram. 



The evaluation was carried out by using a table 
given by Cloud in which he indicates the mean trans- 
mission under various conditions for spectral intervals 
1 -micron wide. Series of values were grouped together 
that have, for a given spectral interval, various 
amounts of C0 2 in the path and also various amounts 
of nitrogen, and the reduced paths were determined 
from (62). This quantity, here briefly called u as 
explained above, was plotted as abscissa — or, rather, 
log u was so plotted. The points then define a trans- 
missivity curve, t(u). If this is done for several inter- 
vals, a number of such curves which run parallel to 
each other is obtained. If the transmissivity curve is 
indeed the same for different intervals, one can bring 
all these curves to coincidence by a shift along the 
log u scale, each shift defining a mean generalized 
absorption coefficient for the interval considered rela- 
tive to some reference interval. The final result, repre- 
senting a standard transmissivity curve for carbon 
dioxide is shown in fig. 10. The measured points fall so 
nearly upon the adopted curve, within the inevitable 
errors of observation and processing, that the points 
shown were not distinguished according to interval, 
C0 2 path, and N 2 pressure. The process also furnishes 
a curve of the generalized absorption coefficient to 
which we shall return. 
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Fig. 10. Transmissivity curve for carbon dioxide. 
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Some further smoothing of the transmissivity curve 
was done numerically when the carbon dioxide radia- 
tion chart was constructed. This served to regularize 
first and second differences. It is important to do this 
smoothing, since later on flux-divergence charts are to 
be obtained by numerical differentiation with respect 
to u\ these become unnecessarily irregular unless the 
curves to be differentiated are smooth to begin with. 
Obviously, the increase in accuracy is purely formal 
and not physical. The flux transmissivity t f was com- 
puted numerically by using (45). For use in a flux- 
divergence chart, the logarithmic derivative, dr/d 
log u, was computed numerically, and the same was 
done for the flux transmissivity, tf. The two trans- 
missivity curves and their logarithmic derivatives are 
shown in fig. 11. It will be seen that r and tf run very 
closely parallel, as implied by the approximation (46), 
and the same is true for the corresponding derivatives. 
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Fig. 12. Generalized absorption coefficient for 15-m 
carbon-dioxide band. 



We found later that the flux-divergence chart can be 
constructed more efficiently and reliably by differentia- 
ting only the completed radiation chart with respect to 
u. Hence, the function drp/d log u was not actually 
used later on in the calculations. 

The generalized absorption coefficient derived from 
the shift of the r-curves for various spectral intervals 
is shown in fig. 12. The main curve, drawn somewhat 
heavier, is the one obtained from the laboratory 
measurements ; it is assumed to correspond to a tem- 
perature of 20C. The lower part of the curve, however, 
on both edges of the band, is a rather arbitrary extra- 
polation. This curve of L was now approximated on 
each side of the band center (667 cm" 1 , from spectro- 
scopic data) by a half -parabola, using a trial-and -error 
procedure. The coefficients of equation (59) are 
b = 0.10; for the short-wave side a = 3.4 • 10~ 4 , and 
for the long-wave side a = 4.6 • 10~ 4 . From this, the 
corresponding band contours for other temperatures 
were computed by using formula (61). It may be noted 
that the point b = 0.10, representing the apex of the 
two half-parabolas, is considerably lower than the 
apex of the actual L-curve for 20 deg which is at 0.44 
in fig. 12. This discrepancy is due to the fact that the 
representative parabolas were chosen so as to match 
the L-curve over a reasonable fraction of the band 
contour rather than just at the apex. 

Fig. 13 shows the actual absorption, 1 — r(uL), 
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FlG. 13. Actual absorption of 15-m carbon-dioxide band. Corresponds to curve marked "trial" in fig. 14. 

curves of fig. 13 are planimetered and the resulting 
points joined, one obtains the curve marked "trial 
curve." It is somewhat too high for long paths (right 
upper corner of fig. 14) but appears too low compared 
to the Ohio State measurements at short paths. When 
our work was almost finished, a paper by Yamamoto 
and Sasamori appeared (1958) in which, for short 
paths, measurements by Kaplan and Eggers (1956) 
are also taken into account. The final curve of band 
areas, then, is likely to lie somewhere near the Ohio 
State values and near Yamamoto's curve. The curve 
ultimately adopted, as shown in fig. 14, was deliber- 
ately taken somewhat low for the following reason. 
Short paths correspond usually to data taken high up 
in the atmosphere ; that is, we want to get our maxi- 
mum accuracy when computing cooling of the strato- 
sphere. Now the experimental data of Howard, Burch, 
and Williams indicate that the pressure dependence 
of the absorption is not quite linear in p but somewhat 
slower, approximately as £ 0 - 88 . Thus, a linear pressure 
correction as applied in practice tends to underestimate 
the band area, and we have tried to compensate for 



computed from the values of r given in fig. 10 and 
those of log L for 20C given in fig. 12. Successive curves 
differ by a factor of ten in u. To the degree in which 
our method of approximation is accurate, these curves 
should duplicate the actual observations of the type 
shown in fig. 9. We next made a further check of figs. 
10 and 12 against the observations, by using the data 
obtained independently by Howard, Burch, and 
Williams at Ohio State University (1955). The latter 
authors used an absorption cell 22 m long in which by 
means of mirrors multiple traversals could be achieved. 
For measurements in the 15-/x carbon-dioxide band, 
4, 8, and 16 traversals were used. Nitrogen was used 
as a broadening agent. Howard, Birch, and Williams 
give their results in the form of band areas. This is the 
area under a given absorptivity curve (maximum 
absorptivity = 1) with wave number in cm -1 as 
abscissa. In order to make the measurements of Cloud 
comparable with these results, we planimetered his 
curves and converted from the wavelength scale to a 
wave-number scale. The results of both sets of observa- 
tions are shown in fig. 14. If now the areas under the 
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FlG. 14. Computed and measured band areas for carbon dioxide, 15-/u band. 



this by shifting the adopted curve to the right in the 
lower-left corner of fig. 14. 

The computation of the C0 2 radiation chart as 
described in the next section was based on the general- 
ized absorption coefficient of fig. 12 and hence corre- 
sponds to the dashed curve of fig. 14. The chart needs 
now to be corrected so as to correspond to the solid 
curve. A simple approximation is suggested by the fact 
that a band covers only a comparatively narrow region 
of the spectrum; hence, the black-body intensity 
should not change very appreciably by function of 
frequency over the band. If the black-body intensity 
was independent of v over the band, the ordinates of 
a radiation chart for given u would be directly propor- 
tional to the band area. Hence, we corrected the 
finished radiation chart in the following way : for any 
value of u (and all temperatures involved), we multi- 
plied the chart ordinates by A/ A i where A is the band 
area (abscissa of fig. 14) corresponding to the adopted 
solid curve and A x is the band area corresponding to 
the dashed trial curve. Since, however, fig. 14 refers 
to beam transmission and the chart to slab trans- 
mission, the correction for a given log u in the chart 
is read off fig. 14 at the argument value log u + 0.20, 
according to equation (46). 



It should be noted in particular that the absorption 
contours of fig. 13 have not been corrected; they do 
correspond to the dashed curve of fig. 14. If necessary, 
they may also be corrected by relabeling the curves of 
fig. 13 with different values of u. This can readily be 
done by using again the two curves of fig. 14. 

There exist quite a number of other bands of CO 2 
in the infrared, but their importance for atmospheric 
radiation transfer is small if not entirely negligible. A 
very strong band is centered at 2350 cm" 1 or a wave- 
length of 4.26 Ae. It represents the fundamental of the 
normal vibration *> 2 (see section III) whereas the 15-m 
band is the fundamental of the normal mode vz. This 
band is situated in a region of the spectrum where for 
normal atmospheric temperatures the black-body 
intensity is very small. Only in very hot air such as 
may be found at mid-day in tropical climates does the 
black-body curve extend appreciably into the region 
covered by this band. In such cases, radiative transfer 
caused by this band may become of some importance. 
Such conditions are, however, sufficiently rare so that 
in all probability these effects should be very small in 
climatological averages. We have therefore omitted a 
separate treatment of this band. 

In addition to the two strong bands representing 
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fundamentals of the C0 2 molecule, there is a large 
number of much weaker bands corresponding to 
higher harmonics of the normal modes. All these may 
be found tabulated in Herzberg's book. The only ones 
that should be of importance for atmospheric radiation 
problems are those that begin (or terminate, as the 
case may be) in the vibrational ground state (0, 0, 0) 
of the molecule, since in all bands where the lower 
vibrational state differs from the ground state there 
will be no appreciable population of the lower state 
at temperatures prevailing in the atmosphere. This 
criterion eliminates all the weak bands on the long- 
wave side of 4 ju. It appears, therefore, that these 
minor bands do not play a significant role in the 
radiative transfer of the atmosphere. 

3. Ozone 

The procedure followed to construct a radiation 
chart for ozone is perfectly analogous in all details to 
the one used for carbon dioxide; hence, we can be 
short in our description. The laboratory transmis- 
sivities used to begin with were those of M. Summer- 
field (1941). The absorption cell was 90 cm long; the which is known to be pressure independent. Summer- 
length of the ozone column reduced to STP ranged field's data were used by Elsasser and King (1953) as 
from 0.2 mm to 5 mm, the latter being considerably in the case of C0 2 to obtain the transmissivity curve 
in excess of the amount usually found in a vertical for reduced path which is shown in fig. 15. This differs 
path in the atmosphere. Oxygen was used as a broaden- clearly from the corresponding curve for carbon diox- 
ing agent. The amount of ozone in the path was deter- ide, fig. 10, not only in shape but also in that there is 
mined by means of the continuous absorption in the a much larger scatter of the observed points, indi- 
Hartley band in the ultraviolet part of the spectrum eating that the linear pressure correction, up/p 0f which 
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Fig. 17. Actual absorption for 9.6-m ozone band. Corresponds to curve marked " computed" in fig. 18. 



was a very good approximation for carbon dioxide is 
a rather poor one for ozone. The curve of fig. 15 was 
subjected again to some numerical smoothing so as 
to regularize first and second differences. Then a flux- 
transmissivity curve, tf, was computed by numerical 
integration by using (45). 

The generalized absorption coefficient derived from 
Summerfield's data is shown in fig. 16. The main band 
is the 9.6-m band which, as indicated in section III, 
may be interpreted as the fundamental of the normal 
mode vz. The exact band center is at 1043 cm" 1 . On 
the short-wave side, there is a weaker companion 
centered at 1110 cm" 1 (9.0 p) which is the fundamental 
pi. In order to obtain the band contour at lower tem- 
peratures by the method described, the following 
values for the parameters of equation (59) were used. 
For the main, 9.6-/* band : b = 0.40, on the short-wave 
side a = 14* 10 -4 , and on the long-wave side a = 
4-10" 4 ; for the satellite band at 9 n\ b = 1, on the 
short-wave side a = 40-10 -4 , and on the long-wave 
side a — 14 -10 -4 . Some of the curves for lower tem- 
perature are shown on fig. 16. 

After the construction of the ozone radiation chart 



had been completed, we became aware of a new and 
very extensive set of measurements by Walshaw(1957) 
which had just been published. Walshaw's absorption 
tube is 35 cm long with a double traversal by the beam. 
The amount of ozone was determined from the absorp- 
tion of an ultraviolet band known to be pressure 
independent. Air was used as the broadening agent of 
the infrared band. Self-broadening of ozone was found 
to be 1.61 times the broadening produced by the same 
amount of air, and this figure was used to reduce all 
data to the condition of infinite dilution of ozone in 
air. Walshaw presents his results in terms of band 
areas under the 9.6-m absorption band. In determining 
these areas, however, he eliminates the 9-jtx band 
altogether by extrapolating the 9.6-ju band graphically 
across the area where the two bands overlap. In order 
to make our results comparable to those of Walshaw, 
we did the same to the actual absorption contours, 
1 — t(uL), for 20C, computed by means of the data 
of fig. 15 and fig. 16. The result, after elimination of 
the 9-m band by graphical means, is shown in fig. 17. 

Fig. 18 shows a plot of Walshaw's data, using our 
standard pressure reduction, plotted against the band 
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Fig. 18. Computed and observed band areas for 9.6-ju ozone band. 



area of the 9.6-^u band. The values derived from the 
band areas of fig. 17 are given by the dashed curve 
marked "computed." Walshaw's measurements were 
made in series of runs, each run containing an approxi- 
mately constant amount of ozone. Points of the same 
run are connected by dashed lines, the figure at the 
top of each line giving an approximate mean value of 
the ozone column in cm at STP. The pressure of the 
broadening air increases upward in each run. The 
points for which the effective air pressure (actual air 
pressure plus equivalent air pressure corresponding to 
self-broadening) is less than 10 cm Hg are marked by 
open circles; the points for which the effective air 
pressure is more than 10 cm Hg are marked by full 
circles. It is at once clear that the linear pressure 
reduction scheme used in our work is a poor approxi- 
mation : the points are very far from clustering about 
a single curve as they should if this pressure correction 
was valid. We did not, however, feel inclined to intro- 
duce a more complicated pressure correction, since 
this would make the use of a radiation chart quite 



cumbersome. Instead, we decided to make a chart that 
is applicable mainly to the stratosphere where most 
of the ozone is concentrated. The "computed" curve 
of fig. 18 was corrected by reducing all band areas by 
10 per cent; this gives the solid curve. Correspond- 
ingly, all values of the computed radiation chart were 
corrected by multiplying them by a factor of 0.90 
throughout. Again, as in the case of C0 2f it should be 
noted that the absorption contours of fig. 17 have not 
been so corrected. The correction brings the calculated 
curve close to Walshaw's measured points in the region 
of low enough pressures and suitable ozone concentra- 
tions to be representative of the stratosphere. We also 
suggest that one might obtain a tolerable approxima- 
tion to the actual pressure dependence on using a 
linear pressure reduction, p/po, up to some fixed 
pressure p u while at higher air pressures using a fixed 
pressure correction, pi/ pa. This would clearly be 
equivalent to assuming that Walshaw 's curves coincide 
with our solid curve for broadening pressures lower 
than p i, whereas for broadening pressures above pi 
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the absorption would be pressure-independent, corres- 
ponding to a vertical course of Walshaw's curves. The 
latter, as fig. 18 shows, is not too far from the truth. 

There is one other absorption band of ozone in the 
far infrared, centered at 710 cm" 1 , or very nearly 14 /x. 
This band is considerably less intense than the 9.6-ju 
band. The analysis of its contribution to radiative 
transfer is complicated by the fact that it is superposed 
upon the very strong 15-/x carbon-dioxide band. In all 
likelihood, however, the band is only of minor impor- 
tance for radiative transfer. For possible future use, 
we give in table 9 a rough estimate of the generalized 
absorption coefficient taken from older measurements 
of Hettner, Pohlmann, and Schumacher (reproduced 
in Herzberg's book). 

4. Water vapor 

The far infrared spectrum ot water vapor is more 
extensive than that of carbon dioxide or ozone. It 



consists of two bands, the 6.3-ja vibration -rotation 
band which is the fundamental of the normal mode z> 2 
and the purely rotational band at the long-wave end of 
the spectrum. The region between the two, in the 
neighborhood of 10 n, is generally known as the 
' 'window' 1 and will be so designated throughout. 

It so happens that the region of the window is the 
least well known part of the water- vapor spectrum, 
whereas the two bands surrounding the window are 
now reasonably well known in respect to their inten- 
sity. We have therefore chosen to compute radiation 
charts for the two bands separately and only at the 
end a radiation chart for the entire water- vapor spec- 
trum. This should facilitate the construction of im- 
proved radiation charts for the overall spectrum when 
further data on the window have become known. 

The first water-vapor radiation chart, which was in 
fact the first radiation nomograph of all, was con- 
structed by Mugge and M oiler (1932). The next such 
chart was developed by Elsasser (1942) based on more 
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Moller (1943), Bruinenberg (1946), Robinson (1950), 
Hales (1951), Yamamoto (1952), and Yamamoto and 
Onishi (1953). The present water-vapor radiation 
chart differs from Yamamoto's only in minor details, 
though it is based almost entirely on experimental 
data obtained later than Yamamoto's work. 

The experiments carried out by Howard, Burch, 
and Williams (1955) furnish transmissivity curves for 
water vapor not only for the 6.3-m band, but the same 
curve applies also for five other bands located in the 
near infrared. There is no reason to question that it 
should also be valid for the rotational band where the 
primary data are inadequate for the construction of 
a separate transmissivity curve. Thus, we have used 
the same transmissivity curve for both the 6.3- y. band 
and the rotational band. Howard, Burch, and Williams 
find that the shape of the transmissivity curve depends 
somewhat on pressure. Their two curves, representing 
a total pressure of 74 cm and 12.5 cm, respectively, 
are plotted in fig. 19, together with curves derived by 
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detailed spectroscopic data. Subsequently, such charts 
were designed by a number of authors — notably, 
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Fig. 22. Actual absorption of 6.3-m water-vapor band. Corresponds to curve marked "computed" in fig. 23. 



Yamamoto (1952) and Daw (1956). All curves have 
been shifted along the abscissa so that they coincide 
at the point of 50 per cent absorption. For our com- 
putations, we adopted the curve of Howard, Burch, 
and Williams corresponding to the lower air pressure 
since we felt that it is more important for the water- 
vapor radiation chart to be accurate at somewhat 
higher levels in the atmosphere than right near the 
ground. The differences between the various curves 
are, in any event, small. 

We first discuss the 6.3-m band. Yamamoto (1952) 
has constructed a curve of the generalized absorption 
coefficient using older experimental data. Daw (1956) 
measured the band in great detail by means of the 
equipment at Johns Hopkins University, the experi- 
mental procedures being somewhat similar to those of 
Cloud for C0 2 as outlined above. He constructed a 
curve of the generalized absorption coefficient com- 
patible with his transmissivity curve of fig. 19. The 
curves of Yamamoto and Daw are compared in fig. 20. 
It appears from this that Daw's curve is somewhat 
high, though it is likely that its shape is correct. A 
final curve obtained by shifting Daw's curve down- 
ward, the top of which is sketched in fig. 20, was 



adopted after a trial-and -error procedure so as to give 
a satisfactory fit with the band areas shown in fig. 23, 
below. 

The curve of the generalized absorption coefficient 
was next corrected for temperature according to the 
method given above in (61). The center of the band is 
at 1595 cm -1 , and each side of the band was corrected 
separately, as for C0 2 . The coefficients of (59) are as 
follows: b = 2.22, for the band's short-wave side a = 
1.75 -10~ 5 , and for the long-wave side a = 2.4- 10~ 5 . 
The resultant curves are shown in fig. 21. We can now 
obtain the band areas for a series of reduced optical 
path and for 20C; these are plotted in fig. 22. Next, 
we compared these results with the experimental data 
obtained by Howard, Burch, and Williams who give 
areas under the band for a variety of optical paths and 
broadening pressures. All the points are shown in fig. 
23 together with the curve that goes through the band 
areas represented in fig. 22. It should be noted that the 
observed points which fall conspicuously below the 
curve represent measurements made on water vapor 
without the admixture of air or with exceedingly small 
amounts of air added. Clearly, such points are not 
characteristic of conditions in the actual atmosphere 
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Fig. 23. Computed and measured band areas for 6.3-/x water-vapor band. 




and may safely be disregarded. A corrected band-area 
curve, fitting the observed points, was then adopted, 
as shown in fig. 23. The subsequent procedure was the 
same as described in the case of C0 2 . First, a radiation 
chart was constructed by means of the preliminary 
curve representing Daw's results as evaluated in fig. 
22. These values were then corrected by multiplying 
for any given value of u by the ratio A/Ai, where A is 
the band area corresponding to the adopted curve and 
A! the band area corresponding to the preliminary 
curve at the point log u + 0.20. Again, it should be 
remembered that fig. 22 corresponds to the preliminary 
curve of fig. 23 and that, if the band areas of fig. 22 
are to be used for practical applications, they should 
be corrected correspondingly. 

Next, we come to the rotational band. The technique 
of measurement becomes exceedingly difficult in the 
very far infrared. Earlier investigators, Yamamoto in 
particular, had therefore relied on direct summing of 
the line intensities computed from quantum me- 
chanics. It is possible, for this band, to compute the 
absolute line intensities, provided only the electro- 
static dipole moment of the water molecule has been 
measured. Yamamoto's band contour obtained in this 
fashion is shown in fig. 24. (Disregard for the moment 
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Fig. 25. Generalized absorption coefficient for rotational 
water- vapor band. 

the region of the window.) It has been shifted along 
the «-axis so as to compensate for the shift of Yama- 
moto's transmissivity curve in fig. 21. There are no 
published recent measurements on this band per- 
formed with a wide enough slit so as to average out 
the individual lines. We were most fortunate to obtain 
by personal communication the results of two inde- 
pendent sets of measurements of this type. Dr. C. 
Harvey Palmer, Jr., 1 in 1958, made a series of wide-slit 
observations with the long absorption tube at Johns 
Hopkins University at a variety of air pressures. The 
vertical lines in fig. 24 are located at the centers of the 
respective spectral intervals and indicate the limits of 
the series of points measured. A few points which 
represent pure water vapor not broadened by atmos- 
pheric air fall outside these intervals and have been 
omitted. Professor E. E. Bell of Ohio State University, 
in 1956, made a series of absorption measurements in 
moist air at atmospheric pressure only. These have 
been reported in an interim technical report to the 
Air Force of which Professor Bell kindly provided a 

1 Personal communication of data obtained at the Labora- 
tory of Astrophysics and Physical Meteorology of Johns Hopkins 
University. 



copy. It is seen that, in the spectral region of v > 200 
cm' 1 , the agreement among the observations of Palmer 
and Bell and between these and Yamamoto's calcu- 
lated curve is good. For waves longer than 200 cm -1 , 
the agreement of Bell's measurements with Yama- 
moto's curve is moderate. We did not feel compelled 
to deviate from Yamamoto's computed curve, how- 
ever, given the extreme delicacy of measurements in 
this spectral region. Moreover, in this region the black- 
body intensity is so weak that its contribution to 
atmospheric transfer is quite small. The band contour 
on the short-wave side, beyond about 500 cm -1 , is an 
arbitrary extrapolation. For lower temperatures, a 
band contour was again determined by means of (61). 
The parameters entering into equation (59) are b — 
2.96 and a = 0.98 -10~ 6 . The resulting curves are 
shown in fig. 25. Finally, fig. 26 shows the computed 
band contours for a variety of reduced pressures at 
20C. 

We next come to the region of the window between 
the two water- vapor bands. Absorption there is weak ; 
it is also known that there are a number of very weak 
lines in the interval. There is evidence, on the other 
hand, that much of the absorption comes from the 
' 'wings" of the very strong lines concentrated at the 
peak of the rotational band as well as of the 6.3 ju 
vibration-rotation band. Returning to (14), we see 
that at large distances from the line center the term 
a 2 in the denominator may be neglected as compared 
to (v — i>o) 2 . The net effect of the wings of all the 
distant lines is then a continuous absorption with 

Sioti 

= (64) 

i tt{v — Vi) L 

where the summation extends over all the distant 
lines. This quantity for the window has been computed 
by L. Kaplan (unpublished) and the results have been 
evaluated and compared with some atmospheric 
observations by Roach and Goody (1958). There are 
various reservations to be made. The observations are 
inadequate to establish beyond doubt the accuracy of 
the computations. Moreover, there are theoretical 
reasons for doubting the extrapolation of the Lorentz 
formula (14) to very large distances from the line 
center. Still, at the present time, the values given in 
the paper of Roach and Goody seem to be about the 
best we have, and so they were utilized in our com- 
putations. The scheme employed is shown on the left- 
hand side of fig. 24. In the first place, the extrapolated 
generalized absorption coefficients of the two main 
bands were assumed to give rise to absorption. Upon 
this was superposed the absorption of a continuum 
assumed to extend from 660 cm" 1 to 1220 cm"" 1 . The 
absorption coefficient, k, for the latter is also shown in 
fig. 24. Now this continuous absorption is expressed 
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Fig. 26. Actual absorption of rotational water-vapor band. 



by a transmission function, r = e~ kx , for a straight 
beam ; the corresponding transmission function for a 
slab may be obtained from (45) in the form 



tf(Jz 



u) = 2 f 



r(tu)dZ/£ z = 2Eu(ku) (65) 



where Eiz is the third exponential integral in the usual 
designation. We have plotted this function in fig. 27 
together with the slab-transmissivity function corre- 
sponding to our adopted transmissivity for the water- 
vapor bands, as shown in fig. 19. The function Eiz has, 
however, been shifted to the left by one decadic cycle 
of the logarithmic u scale — that is, by a factor of ten. 
This was done for the following practical reason: 
mostly, the absorption in the window region is rather 
small, and, by shifting the continuum transmissivity 
in this way, the two curves of fig. 27 are made to run 
close to each other in the critical region of small 
absorption, and the difference in their behavior is 
made apparent. Now, however, it becomes necessary 
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Fig. 27. Comparison of slab transmissivities of water- vapor 
bands (corresponds to beam transmissivity of fig. 19) and 
continuum. 
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to shift the curve of the generalized absorption coeffi- 
cient for the continuum correspondingly — namely, 
downward by one cycle in the logarithmic scale of fig. 
24. It is seen that in the absence of this shift the 
absorption coefficient would indeed run very closely 
to Yamamoto's values, only slightly above the latter. 

The total transmission in the window was computed 
from the formula 

Ttotal — Tcont. * Tband. (66) 

A formula such as (66) is rigorously correct if two 
continuous transmissions are superposed. It should 
therefore represent a fair degree of approximation for 
the case where at least one of the components is a 
continuum. A radiation chart for the window region 
was next computed by using rtotai, and the final radia- 



tion chart for water vappr was obtained by adding up 
the charts for the two bands and the window. The 
window region as defined here (660-1220 cm -1 , with 
the terms of the numerical integration centered at 
multiples of 40 cm -1 — i.e., from 680 cm -1 to 1200 cm" 1 ) 
extends into the two adjacent bands. Hence, on adding 
up the three contributions, a correction term has to 
be subtracted, representing the fact that the band 
edges, above 660 cm -1 and below 1220 cm -1 respect- 
ively, have been counted twice. This should be noted 
by those who wish to use our figures for further studies 
of the spectroscopic details of water-vapor absorption. 

All other bands of the water-vapor spectrum besides 
those treated are in the near infrared, in a region 
where the black-body intensity at atmospheric tem- 
peratures is vanishingly small. 



V- THE TRANSFER EQUATIONS, TABLES AND CHARTS 



1. Integration of the transfer equations 

We now consider atmospheres with an arbitrary 
vertical distribution of radiating gases. These atmos- 
pheres will be assumed horizontally stratified — that is, 
u is independent of horizontal coordinates and a func- 
tion of height only. In place of (38), we now write 

p Wl 

du* = du—\—. (67) 

Often, the square-root dependence on temperature is 
rather slight and might be omitted. In practice, the 
integration of (67) is carried out by a numerical 
approximation, dividing the atmospheres into layers 
of finite thickness whose contributions are added up. 
As before, we propose to omit the asterisk in this 
chapter and designate by u the reduced vertical path as 
defined by (67). If u is measured in the vertical, the 
path of an inclined beam through a layer of thickness 
u is u sec # (fig. 5). This becomes infinite for & = tt/2, 
and, in order to avoid this singularity, it is convenient 
to split the transfer equation into two parts, one 
equation referring to beams having an upward com- 
ponent of travel and another equation referring to 
beams having a downward component. In accordance 
with this splitting, we may assume 0< # < tt/2 so that 
sec # is always positive. We may further assume that 
u is always measured in the same direction — say, 
increasing upward. We let the intensity, /, of flux /, 
be positive in the direction of propagation of each 
beam. For an inclined beam, the equation of transfer 



(12) now becomes 

dl/du = - k sec - J»] (up) 
dl/du = k sec #[/ — lb] (down) 

where /& is again the black-body intensity. 

On integrating (68) for the actual atmosphere and 
a given radiating gas, we need to remember that, in 
any one atmosphere, u(z) and T(z) are given functions 
of height, z. The pressure effect on the radiation has 
been taken into account by (67), and the pressure 
need no longer be considered (except of course when 
we want to use it as independent variable in place of z). 
On eliminating s, we can, for any one radiating gas, 
characterize a specific atmosphere by the functional 
relationship 

u = u(T) or T = T(u). (69) 

We shall write down the integral of the second 
equation (68) ; this differs from that for the first only 
by the algebraic sign. Equations (68) refer to an 
infinitesimal spectral interval. For such an interval, 
now, let 

r(u) = er hum *. (70) 

The integral of the second equation (68) giving the 
downward flux at level Uo is 

C ul dr(u - Uo) 

I(u 0 ) - - /*(«)— ~du (71) 

Juq du 

as may readily be verified by direct substitution by 
using (70). Ib is a function of temperature only, and 
I b (u) is defined by (69). At the upper boundary, Ui = 
0, we have 1(0) = 0. Sometimes we are concerned 
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with an atmospheric layer of finite thickness bounded 
from above by a cloud. The cloud surface may be con- 
sidered a black body for all practical purposes. Let us 
assume for the sake of generality that the temperature 
of the cloud bottom is not identical with that of the 
air adjacent to it; we may write I h (bd) for the black- 
body radiation emerging from this boundary. If we 
let u = uo at the boundary, the solution fulfilling the 
correct boundary condition is 



I(uo) = I b (bd)r(ui - uo) 



/ 



1,1 dr(u - Uo) 

h(u)— -du. (72) 

du 



Formula (72) becomes a solution of the first equation 
(68), applying to an upward beam, by changing the 
sign of the integral. It is to be remembered that / is 
always positive in our definition. The lower boundary 
may now be the ground, or else a cloud top, usually 
radiating as a black body. 

For simplicity, let us next confine ourselves to (71) 
where there is no separate boundary term. We can 
perform two linear operations on this equation — 
namely, the integration over the spectrum and the 
integration over the angles. 

The integration over the spectrum is carried out in 
two steps. First, we go from the monochromatic trans- 
missivity (70) to the transmissivity averaged over a 
finite interval Av— that is, to r(uL) where L is the 
generalized absorption coefficient. This r is defined by 
(16) and, on suitable normalization of i" (or else f) per 
unit spectral interval, we may set Av = 1. Under this 
operation, (71) remains formally unchanged. There- 
after, we can carry out the full integration over the 
spectrum. 

By integrating next over the polar angles, <p f we 
found in section 1 1. 6 that this does not change the 
basic transfer formalism; we merely have to replace 
the intensity J by the flux / and replace r by the flux 
transmissivity tf- We write B v for the black-body flux, 
where by (42) we have B v = irl h . 

Now (71) becomes, if F designates the integrated 
flux, 



F(u 0 ) = - f dv f 1 5, 

J Juo 



Ul 3tf(u — uo) 



du 



du. (73) 



A similar transformation applies to (72) and the 
corresponding formula for the upward flux. Equation 
(73) and similar ones with boundary terms obtained 
by integrating (72) over the spectrum and angles 
solve the radiative-transfer problem. 

2. Construction of tables and charts 

Equation (73) may be used for the numerical com- 
putation of radiative fluxes. Instead, the process may 



be carried out graphically on a nomogram, a radiation 
chart. Thus, Moller's (1943) radiation chart has log u 
as abscissa and the quantity, say 

P(u,T) = - f Bv (T) d -^d v , (74) 
J du 

as ordinate. Isotherms, curves of constant T, are 
permanently drawn on the chart. By using these 
curves, one can plot the relationship T(u) representing 
a given atmosphere. The flux at u 0 then becomes 

F(u 0 ) = f P(u - uo, T)du. (75) 

Juo 

and may be interpreted as the area under the T(u) 
curve which area is most readily evaluated by means 
of a planimeter. Thus, radiative fluxes can be deter- 
mined rapidly. 

There is clearly a great variety of ways in which the 
basic integrals of the transfer equations can be ob- 
tained numerically or graphically. There seems to be 
little point here in reproducing some of the contents 
of a text on numerical analysis or of one on nomo- 
graphic methods. The transformations between the 
various methods are as a rule straightforward, and 
various patterns will occur to anyone who is confronted 
with the task of evaluating a considerable number of 
radiative fluxes. We shall therefore omit altogether 
the description of numerical techniques of evaluating 
these fluxes and shall instead exemplify the whole 
method in terms of a graphical procedure only. This 
does not mean that this method is superior ; quite often 
the numerical approach might be more convenient, 
but the graphical method is more easily demonstrated. 

In the charts of the author and of Yamamoto, 
among others, a change of variables has been per- 
formed. Instead of u y the temperature, T, is used as 
the independent variable by Elsasser (1942) and the 
black-body flux B by Yamamoto (1952). Subject the 
second integral in (73) to an integration by parts, 

dr F (u) 
Juo du 



B V TF 
B V T F 



* l dB v 

~ I —TF{u)du 
uo Juq CLl 

Tl _ r Tl dB v 

To JT 



ho dT 



r F (u)dT, (76) 



where it is to be remembered that in a given atmos- 
phere the functional relation (69) holds. In analogy 
to (74), we introduce a function of /^variables, 



(77) 
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Fig. 28. Radiation chart. 

integrated over the spectrum. By carrying (76) and 
(77) into (73), we get 

l} T F (u-Uo)dv\ + I Q(u-uq, T)dT, 

To J To 

(78) 



F(w 0 ) = -|J B v7 



where the path of the integral is along the actual 
relationship, u(T) t valid for the atmosphere con- 
sidered. A glance at (77) will show that the two 
boundary terms ahead of the integral in (78) are 
equivalent to integrals of the form 



J 0 



QdT taken along u = const. 



Furthermore, by replacing u — uo by u in the last inte- 
gral in (78), which can be done without loss of gener- 
ality as will presently appear, this formula becomes 



F = f ° 0(0, T)dT + f 1 Q(u(T), T)dT 

+ f°Q(u u T)dT. (79) 

This form shows the integral extended along a closed 
path in the T-Q plane. The integration goes from the 
absolute zero of temperature along the path u — m to 
the temperature To, thence along the actual u(T) 
relationship in the atmospheric layer considered to 7\, 
and from there back to the absolute zero along u = U\. 
The replacement of u — Uo by u in the second integral 
means simply that the zero of optical thickness must, 
in this system of evaluation, be taken at the level for 
which the flux is computed. 

Let us remark that the quantity tabulated in 
section Vi is not Q itself but is 



dB r d3 v 

R(u > T)=—-Q(u t T) = J —{\- T {u))dv. (80) 

The advantage is that R vanishes for all points of the 



spectrum where the substance considered is trans- 
parent, whereas Q ~ dB/dT in this part of the 
spectrum. 

3. Use of a radiation chart 

We shall exemplify the determination of radiative 
fluxes by means of the type of chart used earlier by 
the author (1942). If we take the temperature as 
abscissa and Q as ordinate, we get a diagram in which 
(?max = 4(tZ^ is the upper edge. This is inconvenient 
and we transform to new variables, say 



x = aT 2 



y = Q/2aT. 



(81) 



This transformation preserves areas, and the upper 
edge is now a straight line, giving the chart a triangular 
shape (fig. 28). On the chart are drawn curves of con- 
stant u } isopleths of u as we shall call them, implicitly 
defined by (77). The abscissas are labelled in terms of 
F, increasing from right to left, the ordinates are 
labelled in terms of the u isopleths, and a standard 
area is provided for calibration of the planimeter. We 
may remark here that Yamamoto's chart (1952) is 
very similar but rectangular, the abscissa being propor- 
tional to B = oT 4 and the ordinate to Q/T*. 

A radiative flux is equivalent to an area on the radia- 
tion chart; we shall illustrate this by a series of 
examples of increasing complexity. We shall first assume 
a substance such as water vapor which is opaque over 
the whole of the far infrared spectrum for a layer of 
sufficient thickness. The isopleth u = 0 coincides with 
the upper edge of the chart, whereas the isopleth u = 
oo coincides with the horizontal axis. It is readily seen 
that a triangle bounded on the left by a vertical line 
T — const, is equal to the black-body flux, B = aT 4 . 
This graphical representation embodies the rather 
obvious theorem that an isothermal layer of infinite 
thickness of a substance which is nowhere in the spec- 
trum completely transparent radiates as a black body. 
Next, take the nearly triangular area shaded in fig. 28 
which is bounded from above by an isopleth, u = 
const. It represents the quantity 



f QdT = f 



B v T F (u)dv, 



(82) 



that is, the radiation transmitted by a slab of thickness 
u (itself not radiating) upon which the black body 
flux B impinges. The area above the shaded region 
and bounded on the left by the same vertical line is 
then, 



f RdT = f B t (l - r F {u))dv, 



(83) 



that is, the emission of an isothermal layer of thick- 
ness u and temperature T. All these results can readily 
be derived from (77). 
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Fig. 29. Computation of flux on radiation chart. 



Next consider the emission of a layer with arbitrary 
stratification. This is given by (79), and one readily 
verifies that this flux is represented by the area (a) in 
fig. 29. We recall that the zero of optical thickness 
must here be taken at the level for which the flux is 
determined. At any level in the atmosphere, there will 
be a flux in either direction, and the net flux, counted 
positive if going upward, will be the difference, 

-Fnet = F nj> — Fdown- (84) 

For instance, if the surface considered is at the ground, 
the up-going flux will be the black-body radiation of 
the ground whose temperature may now be TV It 
follows at once that area (b) in fig. 29 represents the 
net flux leaving the ground. Again, assume that we 
want to determine the net flux at a higher level in the 
atmosphere ; now, let To stand for the temperature at 
that level. The significance of area (a) with respect to 
the downward flux is the same as before. We now plot 
the actual relation u(T) going downward from the 
reference level, with the reference level as zero. We 
can proceed in this way until we have reached the 
ground. We may next replace the ground by an iso- 
thermal layer of infinite optical thickness at tempera- 
ture r 2 . This is graphically represented by going down- 
wards along the vertical T 2 to the axis of abscissae. 
This is shown in fig. 29, the net upward flux at the 
reference level of temperature To being now given by 
the area (a) + (b) + (c). Since (a) is the downgoing 
flux, it follows that the net flux (84) at the reference 
level is represented by (b) + (c). In a slightly more 
complicated case, the temperature of the ground, say 
r 3 , may be somewhat different from the temperature 
r 2 of the air immediately adjacent to the ground. This 
corresponds to the earlier equation (72) now integrated 
over the spectrum and the angles. One now goes along 
the actual atmospheric curve u(T) to T 2l then along 




Fig. 30. Computation of cooling of a slab on radiation chart. 

the w-isopleth to T 3 and from there on vertically 
downwards. 

Almost all cloud surfaces may be considered black 
bodies and may be accounted for in the radiation chart 
on replacing them by an isothermal layer of infinite 
thickness. If, for instance, in the case shown in fig. 29, 
there is a cloud base at a level of temperature T\ and 
if the reference level To is between this and the ground, 
the net radiation at T 0 would be (b') + (c) where (b') 
is the truncated part of area (b) to the left of a vertical 
line at 7Y 

The u(T) curve as drawn in fig. 29 assumes that the 
temperature decreases with increasing height, which 
is the usual case in the atmosphere. If temperature 
inversions are present, the curves to be drawn on the 
radiation chart become a little more complicated, and 
closer attention may have to be paid to the algebraic 
sign of the different areas ; however, there is no partic- 
ular difficulty, and we may leave out details that can 
readily be supplied as needed. 

An area such as that on the right of the vertical 
T = Ti in fig. 29 must be evaluated in all radiative 
flux calculations. To facilitate this repetitive process, 
the areas, given by (83) have been evaluated once and 
for all for a sufficiently low temperature. Table 20 
gives the areas for a temperature of — 80C for all 
the bands considered. Having these, it is possible to 
design a truncated radiation chart, leaving out the 
triangular area to the right of — 80C. This is shown 
schematically in fig. 30. The area cross-hatched there 
is the difference of two such areas. 

Let us finally lift the restriction that for sufficient 
thickness the substance in question is completely 
opaque in all parts of the spectrum. This is true for 
water vapor, and the earlier radiation charts applied 
to water vapor. Carbon dioxide and ozone, on the 
other hand, cover only a small part of the far infrared 
spectrum. In this case, it is convenient to define a 
spectral interval within which absorption is every- 
where appreciable at the largest thicknesses encoun- 
tered in the atmosphere. One may then repeat all the 



34 



METEOROLOGICAL MONOGRAPHS 



Vol. 4, No. 23 



formal arguments of this section for the limited inter- 
val envisaged, and one may construct a radiation chart 
whose edge represents the fraction of the black-body 
flux which is contained in the interval selected. For 
the case of carbon-dioxide, we shall revert to this 
question when we deal with the spectral overlap of 
H 2 0 and C0 2 . For ozone in the stratosphere, no prob- 
lem arises once the limiting interval is chosen (e.g., 
970-1130 cm" 1 ), since water vapor may be considered 
transparent in the window region except for the emis- 
sion of layers of large moisture near the ground which, 
on combining it with the radiation of the ground, may, 
relative to the ozone band, be considered as external 
black-body radiation of a suitable representative 
temperature. 

4. Radiative cooling 

Assume now two levels in the atmosphere, say (1) 
and (2), where we shall take (2) to be the higher level. 
With the convention (84) which makes the net flux 
positive upwards, the net loss of radiant energy per 
unit time suffered by the layer between (1) and (2) is 



AF = F n J» - R^Q) 



(85) 



where AF is positive for cooling. The evaluation of 
such an expression on the radiation chart is shown in 
fig. 30, the peak corresponding to level (1) being on 
the left, that to level (2) on the right. The rate of 
cooling averaged over the layer is now 



dT 

pc p Az — = AF 
dt 



(86) 



where p is the mean air density in the layer, c p specific 
heat at constant pressure, and Az the thickness of the 
layer. From the barometric equation, this can be 
written 

dT g AF 

— = -— (87) 
dt c P Ap 

where Ap is the pressure difference between the boun- 
daries of the layer. If AF is in gcal per square cm day, 
if Ap is in millibars, and if / is in days, the numerical 
value of the constant on the right of (87) is 4.09. 

In order to compute radiative cooling in the free 
atmosphere, it is usually adequate to subdivide the 
atmosphere into layers of finite thickness and to use 
such formulas as (87). The more rigorous procedure 
would of course be to go to the limit of infinitesimal 
vertical differences, giving 



dT 
dt 



g_dF 
c p dp 



--( 

c» \ 



dF dT dF du 
dT dp + du dp 



) 



(88) 



However, dF/dT vanishes, as may most readily be seen 
by noting that F = fRdT where the integration with 



respect to T is over a closed contour. Hence, (88) 
becomes 

dT g dudF g dlogu dF 
dt c p dp du 



c p dp d log u 
g dlogu 



dp 



f 



dR 



d log u 



dT (89) 



where the integration is again over a closed contour. 
Thus, the determination of radiative cooling can be 
carried out quite directly by means of a nomogram in 
which the quantity di?/dlogw appears as ordinate. 
This quantity is tabulated in section VI for all three 
gases. These tables are less accurate than those for R 
because the numerical differentiation engenders a 
certain loss of accuracy. The procedure of using these 
''flux divergence charts" is analogous to that for the 
radiation chart outlined above. Such charts have been 
designed for water vapor by Bruinenberg (1946), 
Hales (1951), and Yamamoto and Onishi (1953) 
among others. 

5. The carbon-dioxide — water overlap 

In the region of the spectrum where carbon dioxide 
absorbs, there is also some water- vapor absorption. 
A detailed quantitative treatment of this effect upon 
atmospheric radiation would be extremely laborious, 
and we must confine ourselves to a rather crude 
method. This is not too serious, because the spectral 
region involved here contributes at most about 18 per 
cent of the total black-body flux, and often much less. 
Thus, a moderate percentage error in the overlap effect 
corresponds to a comparatively small error in the 
total fluxes. We shall follow the method used by 
Yamamoto (1952). 

The approximation is essentially based upon the 
assumption, already used in (66), that the total trans- 
missivity is the product of the transmissivities corre- 
sponding to the two gases. In other words, in the 
interval of overlap, we replace tf by tf(H20)tf(C02). 
This product rule is of course rigorously correct only 
for the superposition of continuous spectra and is 
inadequate for band spectra. Nevertheless, one can 
expect to get a moderately good approximation owing 
to the fact that most of the absorption takes place in 
the wings of the lines which act somewhat like a con- 
tinuum (the case of absorption proportional to the 
square root of the optical path). We may apply this 
correction either to a water-vapor calculation or else 
to a carbon-dioxide calculation. Consider first the 
former. Referring to the definition (77) of Q, we 
replace it by 



Qcorr = Q - AQ, 



Reorr ~ R — AR 



(90) 
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where now 

fdB v 
tf(H 2 0)[1 - T F {CO*)1dv (91) 
dT 

which, with the integration taken over the spectral 
interval of overlap, effects the replacement considered 
in this interval. Similarly, of course, if we want to 
correct CO2 radiation for underlying water-vapor 
absorption, we write 

AQ - - AR = f — r F (C0 2 )[l - T F (U,0)2dp. (92) 
J dT 

In the tables of these quantities given in the appendix, 
we have taken the spectral interval as 540-820 cm" 1 . 
This particular choice of interval resulted from our 
replacing the integrations in (91) and (92) by simple 
summations over intervals 40 cm" 1 wide, with dB v /dT 
taken at the centers of the intervals. The first and the 
last centers were 560 cm" 1 and 800 cm -1 , respectively. 
The numerical computation of (91) or (92) would still 
be prohibitively lengthy since there are three param- 
eters, the temperature and the optical paths of both 
water vapor and carbon dioxide. By a fortunate cir- 
cumstance, the proportion of the black-body radiation 
contained in this interval relative to the total black- 
body radiation changes rather little with temperature 
(except at the very lowest temperatures, since the 
fraction must clearly go to zero for T ~ 0). Since the 
overlap is in practice most important in the lowest 
layers of the atmosphere, the calculations have been 
carried through for a temperature of 20C throughout. 
The values of (91) and (92) are given in tables 21 and 
22. 

The computation, numerical or graphical, of the 
correction term f AQdt proceeds in exactly the same 
way as the computation of the flux fQdT. Yamamoto 
presents a graph on which AQ is read off, the abscissae 
and ordinates being the carbon dioxide and water 
vapor paths, respectively. He then has lines AQ — 
const, printed on his water vapor radiation chart; 
these lines are at the bottom of the chart, the horizon- 
tal axis being AQ =0, and the higher values of AQ 
cover a strip above the axis. One then draws a curve 
corresponding to the actual atmosphere, and the net 
flux of water vapor radiation corrected for C0 2 is the 
area between the upper and the lower curves (fig. 31). 

In actual calculations of radiative transfer the appli- 
cation of the Yamamoto corrections (91) and (92) 
might offer some difficulties. Note that this correction 
amounts in (91) essentially to replacing J*B' F(RzO)dv 
by the term % fB f r F (H20)TF(C0 2 )dp within the region 
of the overlap, and conversely in (92). Thus if we 
corrected both the water vapor chart and the carbon 
dioxide chart, we would introduce the cross term (the 
last integral) twice. Notice, however, that water vapor 




Fig. 31. Carbon-dioxide overlap correction arranged at 
bottom of water-vapor radiation chart. 



has a sufficient concentration to absorb or emit in the 
region of the overlap only in the lowest layers, say the 
lowest 1-3 km, of the atmosphere. Thus for the com- 
putation of outgoing radiation and atmospheric cooling 
it might commonly suffice to correct only the C0 2 
chart, leaving the H 2 0 chart uncorrected. In any 
event, if the usual procedure is applied of obtaining 
the black-body radiation of the ground by introducing 
a layer of infinite optical thickness, one must be sure 
that the region covered by the C0 2 band is not counted 
twice in the black-body contribution. Another way of 
applying the overlap correction would consist in taking 
the arithmetical mean of (91) and (92) and applying 
this correction once to the combined result of the H 2 0 
and C0 2 charts. The elaboration of a suitable pro- 
cedure that takes account of the problems posed by 
the coexistence, in one region, of bands originating 
from two different gases, is best done by those who 
are actively engaged in practical calculations of this 
type. The present author has not done such calcula- 
tions in recent years, and therefore has confined him- 
self to presenting such numerical tables as can be 
obtained from an evaluation of laboratory measure- 
ments made separately on each gas. 

REFERENCES 

Anderson, P. W., 1949: Pressure broadening in the microwave 
and infrared regions. Phys. Rev., 76, 647-661. 

Bell, E. E., 1956: Interim engineering report for period July- 
September 1956. Ohio State Univ., Contract AF 33(616)- 
3312 with Wright Air Development Center. 

Born, M., 1933: Optik. Berlin, J. Springer, 615 pp. 

Bruinenberg, A., 1946: A numerical method for the calculation of 
temperature changes by radiation in the free atmosphere {in 
Dutch). Kon. Nederl. Meteor. Inst., Meded. en Verhand. 
Serie B, Deel 1, No. 1, 57 pp. 

Chandrasekhar, S., 1950: Radiative transfer. Oxford, Clarendon 
Press, 393 pp. 

Cloud, W. H., 1952: The 15 micron band of CO2 broadened by 
nitrogen and helium. Prog. Rept., Johns Hopkins Univ., 
Contract N-onr 248-01, 56 pp. 



36 



METEOROLOGICAL MONOGRAPHS 



Vol. 4, No. 23 



Cohen, E. R., J. W. M. DuMond, T, W. Layton, and J. S. Rollett, 
1955 : Analysis of variance of the 1952 data on the atomic 
constants and a new adjustment, 1955. Revs. Mod. Phys., 
27, 363-380. 

Daw, H. A., 1956: Transmission of radiation through water vapor 
subject to pressure broadening in the region 4.2 microns to 23 
microns. Tech. Rep. No. 10, Univ. Utah, Contract AF 
19(122)-392 with AFCRC, 70 pp. 

Dorsey, N. E., 1940: Properties of ordinary water substance. 
New York, Reinhold Publ. Co., 673 pp (see p. 298). 

Elsasser, W. M., 1938: Mean absorption and equivalent absorp- 
tion coefficient of a band spectrum. Phys. Rev., 54, 126-129. 

Elsasser, W. M., 1942: Heat transfer by infrared radiation in the 
atmosphere. Harvard Meteor. Studies, No. 6, Cambridge, 
Harvard Univ. Press, 107 pp. 

Elsasser, W. M., and J. I. King, 1953: Transmission data for the 
far infrared bands of carbon dioxide and ozone. Tech. Rep. 
No. 9, Univ. Utah, Contract AF 19(122)-392 with AFCRC, 
44 pp. 

Goody, R. M., 1952 : A statistical model for water vapor absorp- 
tion. Quart. J. r. Meteor. Soc, 78, 165-169. 

Goody, R. M., 1954: The physics of the stratosphere. Cambridge, 
England, Cambridge Univ. Press, 87 pp. 

Hales, J. V., 1951: An atmospheric radiation flux divergence chart 
and meridional cross sections of water vapor radiational heat 
losses computed through its use. Thesis, unpubl., Univ. 
California at Los Angeles, 90 pp. 

Herzberg, G., 1945: Molecular spectra and molecular structure, 
Vol. II, infrared and raman spectra. New York, Van 
Nostrand, 632 pp. 

Howard, J. N., D. L. Burch, and D. Williams, 1956: Infrared 
transmission of synthetic atmospheres. /. Opt. Soc. Amer., 
46, 186-189, 237-247, 334-338, 452-455. (Further data 
in: Geophys. Res. Pap. No. 40, AFCRC, November 1955, 
245 pp.) 

Hughes, R. H., 1953: The microwave spectrum and structure of 

ozone. J. Chem. Phys., 21, 959-960. 
Kaplan, L. D., 1950: Line intensities and absorption for the 15 

micron carbon dioxide band. Chem. Phys., 18, 186-189. 
Kaplan, L. D., 1952: On the pressure dependence of radiative 

heat transfer in the atmosphere. J. Meteor., 9, 1-12. 
Kaplan, L. D., 1953: A quasi-static approach to the calculation of 

atmospheric transmission. Proc. Toronto Meteor. Conf., 

43-48. 

Kaplan, L. D., and D. F. Eggers, 1956: Intensity and line width 
of the 15 micron CO2 band determined by a curve-of- 
growth method. J. Chem. Phys., 25, 876-883. 

Kaplan, L. D., M. V. Migeotte, and L. Neven, 1956: 9.6 micron 
band of telluric ozone and its rotational analysis. Chem. 
Phys., 24, 1183-1186. 



Ladenburg, R., and F. Reiche, 1911 : Ueber Selektive Absorption. 

Ann. der Phys., 42, 181-203. 
Mitchell, A. C, and M. K. Zemansky, 1934: Resonance radiation 

and excited atoms. New York, MacMillan, 385 pp. 
Moller, F., 1943: Das strahlungsdiagramm. Berlin, Springer, 9 pp. 
Miigge, R., and F. Moller, 1932 : Zur berechnung von Strahlungs- 

stromen und Temperaturanderungen in Atmospharen von 

beliebigem Aufbau. Z.f. Geophys., 8, 53-64. 
Plass, G. N., 1958: Models for spectral band absorption. /. Opt. 

Soc. Amer., 48, 690-703. 
Plass, G. N., and D. Warner, 1952: Influence of line shift and 

asymmetry of spectral lines on atmospheric heat transfer. 

/. Meteor., 9, 333-339. 
Plass, G. N., and D. I. Fivel, 1953: Influence of Doppler effect and 

damping on line-absorption coefficient and atmospheric 

radiation transfer. Astrophys. J., 117, 225-233. 
Plass, G. N., and D. I. Fivel, 1955: A method for the integration 

of the radiative transfer equation. /. Meteor., 12, 191-200. 
Randall, H. M., D. M. Dennison, N. Ginsburg, and L. R. Weber, 

1937: The far infrared spectrum of water vapor. Phys. 

Rev., 52, 160-174. 
Roach, W. T., and R. M. Goody, 1958: Absorption and emission 

in the atmospheric window from 770 to 1250 cm -1 . Quart. 

J. r. Meteor. Soc, 84, 319-333. 
Robinson, G. D., 1950: Note on the measurement and estimation 

of atmospheric radiation-2. Quart. J. r. Meteor. Soc, 76, 

37-51. 

Strong, J., and G. N. Plass, 1950: The effect of pressure broaden- 
ing of spectral lines on atmospheric temperature. Astro- 
phys. J., 112, 365-379. 

Summerfield, M., 1941 : Pressure dependence of the absorption in 
the 9.6 micron band of ozone. Thesis, Calif. Inst. Tech., 
unpubl., 65 pp. 

Van Vleck, J. H., and V. F. Weisskopf, 1945: On the shape of 

collision-broadened lines. Revs. Mod. Phys., 17, 227-236. 
Walshaw, C. D., 1957: Integrated absorption by the 9.6 ju band 

of ozone. Quart. J. r. Meteor. Soc, 83, 315-321. 
Whittacker, E. T., and G. N. Watson, 1940: Modern analysis, 

fourth ed. Cambridge, England, Cambridge Univ. Press, 

595 pp. 

Wilson, M. K., and R. M. Badger, 1948: A reinvestigation of the 
vibration spectrum of ozone. J. Chem. Phys., 16, 74. 

Yamamoto, G., 1952: On a radiation chart. Sci. Rep. Tdhoku 
Univ., Series 5, Geophys., 4, 9-23. 

Yamamoto, G., and G. Onishi, 1953 : A chart for the calculation of 
radiative temperature changes. Sci. Rep. Tdhoku Univ., 
Series 5, Geophys. 4, 108-115. 

Yamamoto, G., and T. Sasamori, 1958: Calculation of the ab- 
sorption of the 15 ft carbon-dioxide band. Sci. Rep. 
Tdhoku Univ., Series 5, Geophys., 10, 37-45. 



VI. NUMERICAL TABLES 



For definitions, refer to the following formulas : 

dBv/dT (6) computed from Planck formula, 

t (16) determined empirically for each gas, 

tf (45) by numerical integration from r, 

L (39) determined empirically for each 

band, 



R (80) R(u,T) = / ~ r P )dv, 

fRdT. . . (83) fRdT = fB v {\ - r F )dv y and 

(m = const.) (u =const.) 

- AR. (90) (91) overlap corrections H 2 0 - C0 2 . 

The unit of heat flow is gcal per sq cm day through- 
out. The optical path, w, stands for the "reduced" 
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path up/po where po is standard pressure; otherwise, 
use the more accurate reduction formula (67). For 
carbon dioxide and ozone, the vertical thickness u is 
measured by the height in cm of a column reduced to 



standard pressure and temperature. For water vapor, 
it is measured in g per sq cm, or the height of the 
equivalent column of precipitable water. All logarithms 
used in this work are to base ten. 



Table 1. (dBv/dT) -10*. 



v (cm -1 ) 


—80° 


—70° 


—60° 


—50° 


—40° 


—30° 


—20° 


— 10° 


0° 


+ 10° 


+20° 


+30° 


+40° 


40 


85 


85 


85 


85 


85 


85 


85 


85 


85 


85 


85 


85 


85 


80 


333 


334 


335 


335 


336 


337 


337 


338 


338 


338 


339 


339 


339 


120 


723 


727 


731 


735 


738 


741 


743 


745 


747 


749 


750 


752 


753 


160 


1222 


1235 


1247 


1258 


1267 


1275 


1282 


1289 


1295 


1300 


1305 


1309 


1313 


200 


1790 


1821 


1848 


1872 


1893 


1912 


1929 


1944 


1958 


1971 


1982 


1992 


2001 


240 


2384 


2443 


2495 


2541 


2583 


2620 


2653 


2683 


2710 


2734 


2757 


2777 


2796 


280 


2964 


3062 


3150 


3229 


3300 


3363 


3421 


3473 


3520 


3564 


3603 


3639 


3672 


320 


3493 


3642 


3777 


3899 


4010 


4110 


4201 


4284 


4360 


4429 


4492 


4550 


4603 


360 


3945 


4155 


4347 


4522 


4682 


4829 


4962 


5085 


5198 


5301 


5396 


5483 


5564 


400 


4299 


4577 


4834 


5071 


5290 


5492 


5678 


5849 


6007 


6153 


6288 


6413 


6528 


440 


4546 


4896 


5223 


5529 


5813 


6078 


6324 


6552 


6764 


6960 


7143 


7313 


7471 


480 


4684 


5106 


5505 


5883 


6238 


6571 


6883 


7174 


7447 


7702 


7940 


8162 


8370 


520 


4719 


5208 


5678 


6128 


6555 


6959 


7342 


7702 


8042 


8361 


8660 


8942 


9207 


560 


4661 


5211 


5746 


6265 


6763 


7240 


7694 


8126 


8536 


8924 


9290 


9637 


9964 


600 


4523 


5124 


5717 


6299 


6865 


7411 


7937 


8441 


8923 


9383 


9820 


10235 


10629 


640 


4320 


4961 


5603 


6240 


6866 


7478 


8073 


8648 


9201 


9733 


10242 


10729 


11193 


680 


4069 


4737 


5416 


6099 


6778 


7449 


8107 


8748 


9372 


9974 


10555 


11114 


11651 


720 


3784 


4467 


5172 


5889 


6612 


7333 


8047 


8750 


9438 


10109 


10760 


11391 


11999 


760 


3479 


4165 


4883 


5624 


6379 


7141 


7904 


8661 


9408 


10142 


10860 


11560 


12239 


800 


3165 


3843 


4563 


5317 


6094 


6887 


7688 


8491 


9291 


10082 


10861 


11626 


12373 


840 


2852 


3513 


4225 


4980 


5769 


6582 


7412 


8252 


9095 


9937 


10771 


11595 


12406 


880 


2547 


3183 


3879 


4626 


5415 


6239 


7088 


7955 


8833 


9717 


10599 


11476 


12345 


920 


2257 


2861 


3533 


4263 


5044 


5868 


6727 


7612 


8516 


9432 


10355 


11278 


12198 


960 


1985 


2554 


3194 


3900 


4665 


5481 


6339 


7233 


8153 


9094 


10048 


11010 


11974 


1000 


1734 


2264 


2869 


3545 


4286 


5085 


5936 


6829 


7757 


8713 


9690 


10682 


11683 


1040 


1506 


1994 


2560 


3202 


3914 


4690 


5524 


6409 


7336 


8299 


9291 


10304 


11334 


1080 


1299 


1746 


2272 


2875 


3553 


4301 


5112 


5981 


6900 


7862 


8860 


9887 


10936 


1120 


1115 


1521 


2005 


2568 


3209 


3923 


4706 


5552 


6455 


7409 


8405 


9438 


10500 


1160 


953 


1318 


1761 


2283 


2883 


3561 


4311 


5129 


6011 


6948 


7936 


8967 


10034 


1200 


810 


1137 


1539 


2019 


2579 


3217 


3931 


4717 


5571 


6487 


7459 


8481 


9547 


1240 


685 


976 


1339 


1779 


2296 


2894 


3569 


4319 


5142 


6031 


6982 


7989 


9045 


1280 


577 


835 


1161 


1560 


2037 


2592 


3227 


3939 


4726 


5584 


6509 


7495 


8536 


1320 


485 


711 


1002 


1363 


1799 


2313 


2906 


3578 


4328 


5151 


6045 


7005 


8025 


1360 


405 


603 


862 


1187 


1584 


2057 


2608 


3239 


3949 


4735 


5594 


6524 


7519 


1400 


338 


510 


738 


1029 


1389 


1823 


2333 


2922 


3590 


4337 


5160 


6056 


7021 


1440 


281 


430 


631 


890 


1215 


1610 


2080 


2628 


3254 


3960 


4744 


5604 


6536 


1480 


232 


362 


537 


767 


1059 


1417 


1849 


2356 


2941 


3605 


4349 


5170 


6066 


1520 


192 


303 


456 


660 


920 


1244 


1638 


2106 


2650 


3273 


3975 


4755 


5613 


1560 


158 


253 


386 


565 


797 


1090 


1448 


1877 


2381 


2963 


3623 


4363 


5181 


1600 


130 


211 


326 


483 


689 


951 


1276 


1669 


2134 


2675 


3294 


3992 


4769 


1640 


106 


175 


275 


412 


594 


829 


1122 


1480 


1908 


2409 


2987 


3644 


4379 


1680 


87 


145 


231 


351 


511 


720 


984 


1310 


1702 


2165 


2703 


3318 


4012 


1720 


71 


120 


194 


297 


439 


624 


861 


1156 


1514 


1941 


2440 


3015 


3667 


1760 


57 


99 


162 


252 


376 


540 


752 


1018 


1344 


1736 


2198 


2733 


3345 


1800 


46 


82 


135 


213 


321 


466 


656 


895 


1191 


1550 


1976 


2473 


3045 


1840 


38 


67 


113 


180 


274 


402 


570 


785 


1054 


1381 


1773 


2233 


2766 


1880 


30 


55 


94 


151 


233 


346 


495 


688 


930 


1228 


1587 


2013 


2508 


1920 


24 


45 


78 


127 


198 


297 


429 


601 


819 


1090 


1419 


1811 


2271 


1960 


20 


37 


64 


106 


168 


254 


371 


524 


721 


966 


1266 


1626 


2052 


2000 


16 


30 


53 


89 


142 


217 


320 


457 


633 


855 


1128 


1458 


1851 


2040 


12 


24 


44 


74 


120 


186 


276 


397 


555 


755 


1003 


1305 


1667 


2080 


10 


20 


36 


62 


101 


158 


238 


345 


486 


666 


891 


1167 


1499 


2120 


8 


16 


30 


52 


85 


135 


204 


299 


424 


586 


790 


1041 


1346 


2160 


6 


13 


24 


43 


72 


115 


175 


259 


370 


515 


699 


928 


1207 


2200 


5 


10 


20 


36 


60 


97 


150 


224 


323 


453 


618 


826 


1081 


2240 


4 


8 


16 


29 


51 


82 


129 


193 


281 


397 


546 


734 


966 


2280 


3 


7 


13 


24 


42 


70 


110 


167 


244 


348 


482 


652 


863 


2320 


2 


5 


11 


20 


35 


59 


94 


143 


212 


304 


424 


578 


770 


2360 


2 


4 


9 


17 


29 


50 


80 


123 


184 


266 


373 


512 


686 


2400 


1 


3 


7 


14 


25 


42 


68 


106 


159 


232 


328 


453 


610 
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Lo 




e 

— D 


— 4 


— 3 


—2 


—i 


0 


l 


2 


0.0 




99.91 


97.66 


92.18 


81.70 


60. 


27.27 


3.60 


0.14 


0.1 




99.82 


97.27 


91.41 


80.19 


57.62 


23.89 


2.66 


0 


0.2 




99.70 


96.85 


90.59 


78.56 


54.58 


20.66 


1.96 


0.3 




99.55 


96.39 


89.72 


76.80 


51.42 


17.63 


1.47 




0.4 




99.37 


95.90 


88.79 


74.91 


48.15 


14.84 


1.13 




0.5 




99.16 


95.38 


87.80 


72.88 


44.78 


12.30 


0.89 




0.6 




98.92 


94.82 


86.74 


70.71 


41.33 


10.02 


0.71 




0.7 




98.65 


94.22 


85.61 


68.39 


37.82 


8.01 


0.56 




0.8 


100.00 


98.35 


93.58 


84.40 


65.92 


34.28 


6.27 


0.42 




0.9 


99.97 


98.02 


92.90 


83.10 


63.30 


30.75 


4.80 


0.28 




1.0 


99.91 


97.66 


92.18 


81.70 


60.53 


27.27 


3.60 


0.14 





Log u 



0.0 
0.1 

0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 
1.0 



Log u 



0.0 
0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 
1.0 



Table 3. tf for carbon dioxide. 



Log u 


-6 


-5 


-4 


-3 


-2 


-l 


0 


l 


2 


0.0 




99.52 


96.62 


90.24 


77.93 


53.83 


20.67 


2.19 




0.1 




99.37 


96.15 


89.35 


76.14 


50.69 


17.73 


1.63 




0.2 




99.18 


95.65 


88.40 


74.22 


47.46 


15.00 


1.23 




0.3 




98.97 


95.12 


87.39 


72.17 


44.13 


12.51 


0.94 




0.4 




98.72 


94.55 


86.30 


69.97 


40.74 


10.28 


0.73 




0.5 




98.45 


93.94 


85.14 


67.64 


37.31 


8.30 


0.57 




0.6 




98.14 


93.29 


83.90 


65.15 


33.87 


6.58 


0.43 




0.7 




97.81 


92.60 


82.57 


62.53 


30.44 


5.12 


0.31 




0.8 


99.74 


97.44 


91.86 


81.13 


59.76 


27.07 


3.91 


0.19 




0.9 


99.65 


97.04 


91.08 


79.59 


56.86 


23.80 


2.94 


0.10 




1.0 


99.52 


96.62 


90.24 


77.93 


53.83 


20.67 


2.19 


0.02 





Table 4. t f or ozone. 



Table 6. r for water vapor. 



-4 



-3 



-2 



Log u 



— 3 



— 2 



-1 



100.0 
99.2 
98.4 
97.6 



97.6 
96.8 
96.0 
95.2 
94.4 
93.6 
92.7 
91.7 
90.6 
89.4 
88.1 



88.1 
86.6 
84.9 
83.0 
80.9 
78.6 
76.1 
73.4 
70.5 
67.3 
63.8 



63.8 
60.0 
55.9 
51.5 
46.8 
41.8 
36.6 
31.2 
25.7 
20.2 
15.1 



15.1 
10.8 
7.4 
4.8 
2.9 
1.5 
0.5 
0 



0.0 
0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 



100.00 
99.90 
99.67 
99.32 
98.87 
98.33 
97.72 



97.05 
96.33 
95.56 
94.75 
93.89 
92.98 
92.00 
90.94 
89.78 
88.51 



1.0 97.05 87.11 



87.11 
85.56 
83.84 
81.94 
79.84 
77.53 
75.00 
72.25 
69.28 
66.10 
62.72 



62.72 
59.15 
55.41 
51.52 
47.50 
43.38 
39.19 
34.97 
30.76 
26.61 
22.58 



22.58 
18.74 
15.18 
11.98 
9,20 
6.87 
5.00 
3.57 
2.50 
1.68 
1.03 



Table 5. tf for ozone. 



Table 7. tf for water vapor. 



-5 



-4 



-2 



-1 



Log u 



— 2 



-1 



98.2 
97.4 
96.6 
95.7 



95.7 
94.9 
94.1 
93.2 
92.2 
91.2 
90.0 
88.8 
87.4 
85.8 
84.1 



84.1 
82.2 
80.0 
77.7 
75.1 
72.3 
69.3 
66.1 
62.5 
58.7 
54.6 



54.6 
50.3 
45.7 
40.8 
35.8 
30.7 
25.7 
20.7 
16.0 
11.9 
8.4 



8.4 
5.7 
3.6 
2.1 
1.1 
0.4 
0.1 



0.0 
0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 
1.0 



99.32 
98.98 
98.54 
98.01 
97.40 
96.73 
96.00 
95.22 



95.22 
94.39 
93.50 
92.55 
91.53 
90.42 
89.21 
87.88 
86.42 
84.81 
83.04 



83.04 
81.09 
78.95 
76.61 
74.06 
71.30 
68.33 
65.15 
61.78 
58.23 
54.52 



54.52 
50.67 
46.71 
42.67 
38.59 
34.50 
30.44 
26.46 
22.61 
18.95 
15.54 



15.54 
12.45 
9.73 
7.41 
5.50 
3.98 
2.81 
1.93 
1.28 
0.80 
0.45 



1.03 
0.49 
0 



l 



0.45 
0.19 
0 
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Table 8. Log L for carbon dioxide. 



Table 9. Log L for ozone. 







Log L 




Log L 






Log L 




Log L 14 ft band 


V 




(+20°) 


V 


(+20°) 




V 


(+20°) 


V 




estim. 


560 




0.03-5 


690 


0 




970 


0.8-3 








570 




0.95-5 


700 


0.66-1 




980 


0.6-2 


610 




U.4-5 


580 




0.70-4 


710 


0.28-1 




990 


0.2-1 


630 




0.7—4 


590 




0.37-3 


720 


0.92-4 




1000 


0.6-1 


650 




0.3-3 


600 




0.94-3 


730 


0.50-2 




1010 


0.95-1 


670 




0.7-3 


610 




0.50-2 


740 


0.12-2 




1020 


0.2 


690 




v.y—S 


620 




0.95-2 


750 


0.73-3 




1030 


0.3 


710 




0.7-3 


630 




0.39-1 


760 


0.33-3 




1037 


0.5 


725 






640 




0.79-1 


770 


0.92-4 




1045 


0.9-1 


740 




0.9-2 


650 




0.16 


780 


0.48-4 




1052 


0.5 


760 




0.6-3 


660 




0.38 


790 


0.07-4 




1060 


0.35 


780 




0.2-3 


670 




0.43 


800 


0.61-5 




1070 


0.3-1 


800 




0.6-4 


680 




0.32 


810 


0.17-5 




1080 


0.75-2 


820 




0.8-5 














1090 


0.7-2 




















1100 


0.8-2 




















1110 


-1 




















1120 


0.7-2 




















1130 


0.5-3 
















Table 10. Absorption coefficients for water vapor. 












5.2-ju band 






Window 






Rotational band 






Log L 




Log L 




Log L 


log k 




Log L 




Log L 


V 


(20°) 


V 


(20°) 


V 


(20°) 


(cont.) 


V 


(20°) 


V 


OCt°\ 

v^u ) 


1000 


0.77-5 


1640 


2.08 


680 


0.58-2 


0.22-2 


(0) 


(L10) 


480 


0.60 


1040 


0.31-4 


1680 


2.29 


720 


0.08-2 


0.12-2 


40 


2.04 


520 


U.2o 


1080 


0.85-4 


1720 


2.17 


760 


0.56-3 


0.02-2 


80 


2.47 


560 


a or i 

U.oo— 1 


1120 


0.38-3 


1760 


1.86 


800 


-3 


0.94-3 


120 


2.67 


600 


U.45-1 


1160 


0.93-3 


1800 


1.38 


840 


0.41-4 


0.89-3 


160 


2.69 


640 


U.Uo— 1 


1200 


0.47-2 


1840 


0.91 


880 


0.80-5 


0.86-3 


200 


2.61 


680 


n to i 
U.oo— I 


1240 


0.07-1 


1880 


0.50 


920 


0.20-5 


0.83-3 


240 


2.45 


720 


U.Uo-z 


1280 


0.74-1 


1920 


0.14 


960 


0.25-5 


0.82-3 


280 


2.23 


760 


0.56-3 


1320 


0.35 


1960 


0.79-1 


1000 


0.77-5 


0.82-3 


320 


1.96 


800 




1360 


0.85 


2000 


0.41-1 


1040 


0.31-4 


0.83-3 


360 


1.65 


840 


0.41-4 


1400 


1.27 


2040 


0.04-1 


1080 


0.85-4 


0.85-3 


400 


1.30 


880 


0.80-5 


1440 


1.64 


2080 


0.60-2 


1120 


0.38-3 


0.91-3 


440 


0.95 






1480 


2.05 


2120 


0.12-2 


1160 


0.93-3 


0.98-3 










1520 


2.46 


2160 


0.58-3 


1200 


0.47-2 


0.07-2 










1560 


2.01 


2200 


0.97-4 
















1590 


1.49 


2240 


0.31-4 
















1600 


1.56 


2280 


0.59-5 

















Table 11. R for carbon dioxide. 



Log u 


-80° 


-70° 


-60° 


-50° 


-40° 


-30° 


-20° 


-10° 


0° 


+10° 


+20° 


+30° 


+40° 


_4 


0.004 


0.005 


0.006 


0.007 


0.007 


0.008 


0.009 


0.009 


0.009 


0.010 


0.011 


0.012 


0.012 


-4+0.3 


0.006 


0.007 


0.008 


0.009 


0.010 


0.011 


0.012 


0.013 


0.013 


0.014 


0.015 


0.016 


0.016 


-4+0.7 


0.009 


0.010 


0.012 


0.013 


0.015 


0.016 


0.018 


0.019 


0.020 


0.021 


0.022 


0.023 


0.024 


-3 


0.012 


0.014 


0.016 


0.018 


0.020 


0.022 


0.024 


0.026 


0.028 


0.030 


0.031 


0.032 


0.034 


-3+0.3 


0.017 


0.020 


0.022 


0.025 


0.028 


0.031 


0.033 


0.036 


0.039 


0.042 


0.044 


0.046 


0.048 


-3+0.7 


0.027 


0.031 


0.035 


0.039 


0.044 


0.048 


0.052 


0.056 


0.060 


0.065 


0.069 


0.073 


0.076 


-2 


0.038 


0.043 


0.049 


0.055 


0.062 


0.068 


0.074 


0.079 


0.085 


0.091 


0.097 


0.103 


0.108 


-2+0.3 


0.053 


0.060 


0.069 


0.078 


0.087 


0.096 


0.105 


0.112 


0.121 


0.129 


0.137 


0.145 


0.153 


-2+0.7 


0.083 


0.095 


0.109 


0.123 


0.138 


0.153 


0.167 


0.180 


0.193 


0.206 


0.219 


0.232 


0.245 


-1 


0.117 


0.135 


0.155 


0.175 


0.196 


0.217 


0.237 


0.257 


0.276 


0.295 


0.314 


0.333 


0.351 


-1+0.3 


0.162 


0.188 


0.216 


0.244 


0.273 


0,302 


0.331 


0.360 


0.388 


0.415 


0.442 


0.469 


0.495 


-1+0.7 


0.233 


0.272 


0.313 


0.355 


0.398 


0.441 


0.485 


0.529 


0.572 


0.614 


0.655 


0.695 


0.734 


0 


0.288 


0.337 


0.389 


0.443 


0.498 


0.553 


0.609 


0.665 


0.721 


0.776 


0.829 


0.881 


0.932 


0.3 


0.342 


0.402 


0.465 


0.531 


0.598 


0.666 


0.734 


0.803 


0.872 


0.940 


1.006 


1.071 


1.135 


0.7 


0.412 


0.487 


0.565 


0.646 


0.729 


0.814 


0.899 


0.985 


1.071 


1.157 


1.241 


1.324 


1.406 


+ 1 


0.463 


0.549 


0.638 


0.730 


0.825 


0.922 


1.020 


1.119 


1.218 


1.317 


1.415 


1.512 


1.608 


1.3 


0.513 


0.609 


0.709 


0.812 


0.919 


1.028 


1.139 


1.251 


1.363 


1.475 


1.587 


1.698 


1.808 


1.7 


0.578 


0.687 


0.801 


0.919 


1.042 


1.167 


1.295 


1.424 


1.553 


1.683 


1.813 


1.942 


2.070 


2 


0.626 


0.745 


0.869 


0.998 


1.133 


1.270 


1.410 


1.551 


1.693 


1.836 


1.979 


2.121 


2.262 


2.3 


0.673 


0.802 


0.936 


1.076 


1.222 


1.371 


1.523 


1.676 


1.831 


1.986 


2.141 


2.296 


2.449 


2.7 


0.735 


0.877 


1.024 


1.178 


1.339 


1.503 


1.671 


1.840 


2.011 


2.181 


2.351 


2.521 


2.689 


3 


0.781 


0.932 


1.089 


1.253 


1.425 


1.600 


1.779 


1.960 


2.142 


2.323 


2.503 


2.682 


2.860 


3.3 


0.826 


0.986 


1.153 


1.327 


1.509 


1.695 


1.885 


2.077 


2.269 


2.460 


2.649 


2.836 


3.022 
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Table 12. — dR/d log u for carbon dioxide. 



Log u 


— 80° 


— 70° 


—60° 


—50° 


—40° 


—30° 


—20° 


J0O 


0° 




4-20 


0 

4-30 


4-40 


— 4 


0 00^ 


0 006 


0 007 


O 008 
U.UUo 


u.uuy 


U.U1U 


U.U1U 


f\ ni 1 

U.UI 1 


ft m 1 
U.UI 1 


ft m 1 
U.UI 1 


0.011 


0.012 


0.012 


_4_j_0.3 


0 006 


0 007 


0 008 


0 00Q 


U.U1U 


n nil 
U.UI 1 


n ni o 
U.Ulz 


n m 2 
U.Ulo 


ft m a 
U.U14 


0.014 


0.015 


0.016 


0.016 


_4_l_0 7 


0 000 


0 010 


0 Ol 9 


U.Ulo 


n m 4 

U.U14 


U.UIO 


n ni q 
U.Ulo 


ft nin 
U.UzU 


ft no 1 
U.Uzl 


ft nn 
U.Uzz 


0.024 


0.025 


0.026 


_3 


0 013 


0 01 S 


0 017 


0 010 

u.ui y 


U.UZ1 


n no 2 

U.UZO 


n noA 
U.Uzo 


ft niQ 
U.Uzo 


ft nift 
U.UoU 


ft n"2i 
U.Uoz 


0.034 


0.036 


0.038 


— 3+0.3 


0 010 


0 09? 


n nic 


O 098 


n o^i 


U.Uo4 


n n27 
U.Uo/ 


ft ftAft 

U.U4U 


ft ftA 7 

U.U4o 


ft ftA < 

U.U40 


0.049 


0.052 


0.055 




0 030 


0 035 


0 04.0 


O 04 5 


n n5n 
u.uou 


n nc? 


n n/;n 
U.U0U 


ft ftAA 

U.U04 


ft ftAft 

U.Uov 


ft ftl A 

U.U74 


0.079 


0.084 


0.089 


— 2 


0 04.? 


0 040 


0 056 


0 064 


n n79 
u.u/ z 


u.u/y 


n noA 
U.UoO 


ft ftfti 

u.uyo 


ft 1 nn 
U.1UU 


ft 1 n*7 
U.1U7 


0.113 


0.120 


0.127 


— 2 -1-0 3 


0 050 


0 060 


0 080 


O 001 


U. 1UZ 


n in 
U.l lo 


U.lz4 


n 1 i/i 
U.lo4 


ft 1 A A 

U.144 


n 1 c ? 
U.153 


0.162 


0.171 


0.180 


— 2-1-0 7 


0 006 


0 113 


0 1 30 


O 148 

V. 14:0 


n 1 ^ 

U. 10o 


U.Ioz 


n 1 no 


U.zlo 


U.z33 


ft '"I A ft 

U.z40 


0.265 


0.281 


0.297 


— 1 


0 1 32 


0 1 S5 

U. 1 JO 


0 1 78 


0 909 
u.zuz 


n 996 
U.ZZO 


U.ZoU 


n 07c 
U.z / o 


ft 2ftft 

U.oUU 


U.3z5 


0.349 


0.373 


0.396 


0.419 


— 1 -4-0 3 


0 164. 


0 103 


0 993 


0 953 


0 984 


U.OlO 


n 2 a ft 
U.o4v 


ft 2C1 

U.ooZ 


ft A 1 A 

U.414 


0.445 


0.476 


0.506 


0.536 


— 1 +0.7 


0 1 89 


0 ?1 7 


0 9^3 


0 900 
u.zyu 


A 398 


U.oo/ 


U.4U/ 


ft A A 7 

U.44/ 


ft A Q7 


n CT7 

U.oz / 


0.566 


0.605 


0.643 


q 


0.182 


0.218 


0 9^6 
u.zou 


0 9Q5 


n 335 


U.O /0 


ft A i Q 

U.41o 


ft A A.ft 

U.40U 


ft cni 
U.oUz 


U.54o 


0.587 


0.629 


0.671 


0.3 


0.178 


0.215 


O 953 


0 909 

U.i7Z 


fi 319 


n 21 a 

U.o /4 


U.410 


ft A CO 

U.4o0 


ft cm 
U.oUz 


0.546 


0.590 


0.634 


0.678 


0 7 


0.172 


0.209 


O 947 


0 985 


H 394 


U.oOo 


n a f\l 
U.4U/ 


ft a cn 
U.4oU 


0.494 


0.539 


0.584 


0.630 


0.676 


+ l" 


o!l68 


o!204 


0.241 


0.279 


0.318 


0.358 


0.400 


0.443 


0.487 


0.532 


0.577 


0.623 


0.670 


1.3 


0.165 


0.199 


0 235 


0.273 


0.312 


0.352 


0.394 


0.437 


0 4.80 


0 ^94. 

U.OZ4: 




fl 61 c. 


n 661 

U.OOl 


1.7 


0.161 


0.194 


0.229 


o!266 


o!305 


0.345 


o!386 


0^428 


0.470 


0.513 


0.557 


0.602 


0.646 


2 


0.158 


0.191 


0.225 


0.261 


0.300 


0.340 


0.380 


0.421 


0.462 


0,504 


0.547 


0.590 


0.632 


2.3 


0.156 


0.188 


0.222 


0.257 


0.295 


0.334 


0.374 


0.414 


0.454 


0.494 


0.535 


0.575 


0.614 


2.7 


0.154 


0.185 


0.218 


0.252 


0.288 


0.326 


0.365 


0.404 


0.442 


0.479 


0.515 


0.550 


0.583 


3 


0.152 


0.183 


0.215 


0.248 


0.283 


0.320 


0.358 


0.395 


0.431 


0.465 


0.497 


0.527 


0.555 


3.3 


0.150 


0.181 


0.212 


0.244 


0.278 


0.313 


0.349 


0.384 


0.417 


0.448 


0.476 


0.501 


0.523 



Table 13. R for ozone. 



Log u 


-80° 


-70° 


-60° 


-50° 


-40° 


-30° 


-20° 


-10° 


0° 


4-10° 


4-20° 


4-30° 


4-40° 


-4 


0.006 


0.008 


0.010 


0.013 


0,015 


0.017 


0.020 


0.023 


0.026 


0.029 


0.033 


0.036 


0.039 


-4+0.3 


0.009 


0.012 


0.015 


0.019 


0.023 


0.027 


0.031 


0.036 


0.041 


0.046 


0.052 


0.057 


0.062 


-4+0.7 


0.015 


0.020 


0.026 


0.032 


0.039 


0.046 


0.054 


0.063 


0.072 


0.081 


0.091 


0.100 


0.109 


-3 


0.022 


0.029 


0.038 


0.047 


0.057 


0.068 


0.080 


0.093 


0.107 


0.121 


0.135 


0.149 


0.163 


-3+0.3 


0.032 


0.042 


0.054 


0.067 


0.082 


0.099 


0.117 


0.136 


0.156 


0.176 


0.197 


0.218 


0.239 


-3+0.7 


0.050 


0.065 


0.083 


0.103 


0.127 


0.154 


0.183 


0.213 


0.244 


0.276 


0.309 


0.343 


0.377 


-2 


0.066 


0.086 


0.110 


0.137 


0.169 


0.205 


0.244 


0.284 


0.326 


0.369 


0.414 


0.460 


0.506 


-2+0.3 


0.084 


0.110 


0.141 


0.176 


0.217 


0.263 


0.313 


0.365 


0.419 


0.475 


0.534 


0.594 


0.655 


-2+0.7 


0.109 


0.144 


0.185 


0.232 


0.286 


0.347 


0.412 


0.480 


0.552 


0.627 


0.706 


0.787 


0.869 


-1 


0.127 


0.169 


0.217 


0.273 


0.337 


0.409 


0.485 


0.565 


0.651 


0.740 


0.834 


0.931 


1.029 


-1+0.3 


0.144 


0.192 


0.248 


0.312 


0.386 


0.468 


0.554 


0.646 


0.745 


0.848 


0.956 


1.068 


1.182 


-1+0.7 


0.164 


0.220 


0.285 


0.359 


0.445 


0.538 


0.637 


0.745 


0.859 


0.979 


1.105 


1.234 


1.366 


0 


0.177 


0.238 


0.309 


0.389 


0.482 


0.582 


0.690 


0.808 


0.932 


1.063 


1.200 


1.340 


1.483 


0.3 


0.188 


0.253 


0.329 


0.414 


0.512 


0.618 


0.733 


0.858 


0.990 


1.130 


1.275 


1.424 


1.576 


0.7 


0.200 


0.269 


0.349 


0.439 


0.542 


0.653 


0.774 


0.905 


1.044 


1.191 


1.343 


1.499 


1.658 


+ 1 


0.207 


0.278 


0.360 


0.452 


0.557 


0.670 


0.793 


0.926 


1.067 


1.216 


1.370 


1.528 


1.689 
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Table 15. R for 6.3-ju water-vapor band. 



Log u 


-80° 


-70° 


-60° 


-50° 


-40° 


-30° 


-20° 


-10° 


0° 


+ 10° 


4-20° 


4-30° 


4-40° 


-5 


0.01 


0.01 


0.02 


0.03 


0.04 


0.06 


0.08 


0.10 


0.12 


0.15 


0.19 


0.23 


0.27 


-5+0.3 


0.01 


0.02 


0.03 


0.04 


0.05 


0.07 


0.09 


0.12 


0.15 


0.19 


0.23 


0.28 


0.34 


-5+0.7 


0.01 


0.02 


0.03 


0.05 


0.07 


0.09 


0.12 


0.16 


0.20 


0.26 


0.31 


0.38 


0.45 


-4 


0.02 


0.03 


0.04 


0.06 


0.09 


0.12 


0.15 


0.20 


0.25 


0.32 


0.39 


0.47 


0.56 


-4+0.3 


0.02 


0.03 


0.05 


0.08 


0.11 


0.14 


0.19 


0.25 


0.32 


0.39 


0.48 


0.58 


0.70 


-4+0.7 


0.03 


0.04 


0.07 


0.10 


0.14 


0.19 


0.25 


0.33 


0.42 


0.52 


0.64 


0.77 


0.92 


-3 


0.03 


0.05 


0.08 


0.12 


0.17 


0.23 


0.31 


0.40 


0.51 


0.63 


0.78 


0.94 


1.13 


-3+0.3 


0.04 


0.06 


0.10 


0.14 


0.20 


0.28 


0.37 


0.48 


0.61 


0.77 


0.94 


1.14 


1.37 


-3+0.7 


0.05 


0.08 


0.12 


0.18 


0.25 


0.35 


0.46 


0.61 


0.77 


0.97 


1.19 


1.45 


1.73 


-2 


0.06 


0.09 


0.14 


0.21 


0.29 


0.40 


0.54 


0.70 


0.90 


1.13 


1.39 


1.69 


2.02 


-2+0.3 


0.06 


0.10 


0.16 


0.23 


0.33 


0.46 


0.61 


0.80 


1.02 


1.28 


1.59 


1.93 


2.32 


-2+0.7 


0.07 


0.12 


0.18 


0.27 


0.38 


0.53 


0.71 


0.93 


1.19 


1.49 


1.85 


2.25 


2.70 


-1 


0.08 


0.13 


0.20 


0.30 


0.42 


0.58 


0.78 


1.02 


1.31 


1.65 


2.03 


2.48 


2.99 


-1+0.3 


0.09 


0.14 


0.22 


0.32 


0.46 


0.64 


0.85 


1.11 


1.43 


1.79 


2.22 


2.71 


3.27 


-1+0,7 


0.10 


0.16 


0.25 


0.36 


0.51 


0.71 


0.95 


1.24 


1.58 


1.99 


2.46 


3.01 


3.63 


0 


0.11 


' 0.17 


0.26 


0.39 


0.55 


0.76 


1.02 


1.33 


1.70 


2.14 


2.65 


3.23 


3.90 


0.3 


0.12 


0.19 


0.28 


0.42 


0.59 


0.82 


1.09 


1.42 


1.82 


2.29 


2.83 


3.46 


4.18 


0.7 


0.13 


0.21 


0.31 


0.46 


0.65 


0.89 


1.19 


1.56 


1.99 


2.50 


3.09 


3.78 


4.56 


+ 1 


0.14 


0.22 


0.33 


0.49 


0.69 


0.95 


1.27 


1.66 


2.12 


2.66 


3.29 


4.02 


4.84 


1.3 


0.15 


0.23 


0.36 


0.52 


0.74 


1.01 


1.35 


1.76 


2.25 


2.83 


3.49 


4.26 


5.13 


1.7 


0.16 


0.25 


0.39 


0.57 


0.80 


1.10 


1.46 


1.91 


2.43 


3.05 


3.76 


4.57 


5.50 



Table 16. R for rotational water- vapor band. 



Log u 


-80° 


-70° 


-60° 


-50° 


-40° 


-30° 


-20° 


-10° 


0° 


4-10° 


4-20° 


4-30° 


4-40° 


-5 


0.04 


0.04 


0.04 


0.04 


0.04 


0.05 


0.05 


0.05 


0.05 


0.05 


0.06 


0.06 


0.06 


-5+0.3 


0.06 


0.06 


0.06 


0.07 


0.07 


0.07 


0.08 


0.08 


0.08 


0.08 


0.09 


0.09 


0.09 


-5+0.7 


0.09 


0.10 


0.10 


0.11 


0.11 


0.12 


0.13 


0.13 


0.14 


0.14 


0.15 


0.15 


0.16 


-4 


0.13 


0.14 


0.15 


0.16 


0.16 


0.17 


0.18 


0.19 


0.20 


0.21 


0.21 


0.22 


0.23 


-4+0.3 


0.18 


0.19 


0.20 


0.22 


0.23 


0.24 


0.25 


0.26 


0.28 


0.29 


0.30 


0.31 


0.33 


-4+0.7 


0.26 


0.28 


0.30 


0.32 


0.34 


0.36 


0.38 


0.40 


0.42 


0.44 


0.46 


0.48 


0.50 


-3 


0.34 


0.37 


0.40 


0.42 


0.45 


0.47 


0.50 


0.53 


0.55 


0.58 


0.61 


0.63 


0.66 


-3+0.3 


0.43 


0.47 


0.50 


0.53 


0.57 


0.60 


0.64 


0.67 


0.71 


0.75 


0.78 


0.82 


0.86 


-3+0.7 


0.56 


0.61 


0.65 


0.70 


0.75 


0.80 


0.85 


0.90 


0.95 


1.00 


1.05 


1.11 


1.16 


-2 


0.66 


0.72 


0.77 


0.83 


0.89 


0.95 


1.01 


1.08 


1.14 


1.21 


1.28 


1.35 


1.42 


-2+0.3 


0.76 


0.83 


0.90 


0.97 


1.04 


1.11 


1.19 


1.27 


1.35 


1.43 


1.51 


1.60 


1.69 


-2+0.7 


0.89 


0.98 


1.06 


1.15 


1.24 


1.34 


1.43 


1.53 


1.64 


1.74 


1.85 


1.96 


2.07 


+ 1 


0.99 


1.09 


1.19 


1.29 


1.40 


1.51 


1.62 


1.74 


1.86 


1.98 


2.11 


2.24 


2.37 


-1+0.3 


1.09 


1.20 


1.31 


1.43 


1.55 


1.68 


1.81 


1.95 


2.09 


2.23 


2.37 


2.52 


2.68 


-1+0.7 


1.22 


1.35 


1.48 


1.62 


1.76 


1.91 


2.07 


2.23 


2.39 


2.56 


2.73 


2.91 


3.10 


0 


1.31 


1.46 


1.61 


1.76 


1.92 


2.09 


2.26 


2.44 


2.62 


2.81 


3.01 


3.21 


3.41 


0.3 


1.41 


1.56 


1.73 


1.90 


2.07 


2.26 


2.45 


2.64 


2.85 


3.06 


3.27 


3.50 


3.72 


0.7 


1.53 


1.70 


1.89 


2.07 


2.27 


2.48 


2.69 


2.92 


3.15 


3.38 


3.63 


3.87 


4.13 


+ 1 


1.61 


1.80 


2.00 


2.20 


2.42 


2.64 


2.87 


3.11 


3.36 


3.62 


3.88 


4.15 


4.42 


1.3 


1.70 


1.90 


2.11 


2.33 


2.56 


2.80 


3.04 


3.30 


3.57 


3.84 


4.12 


4.41 


4.70 


1.7 


1.80 


2.02 


2.25 


2.48 


2.73 


2.99 


3.27 


3.55 


3.84 


4.13 


4.43 


4.74 


5.05 



Table 14. —dR/d log u for ozone. 



Log u 


-80° 


-70° 


-60° 


-50° 


-40° 


-30° 


-20° 


-10° 


0° 


+10° 


+20° 


+30° 


+40° 


-4 


0.01 


0.01 


0.02 


0.02 


0.02 


0.03 


0.03 


0.04 


0.04 


0.05 


0.05 


0.06 


0.06 


-4+0.3 


0.01 


0.02 


0.02 


0.03 


0.03 


0.04 


0.04 


0.05 


0.06 


0.07 


0.08 


0.08 


0.09 


-4+0.7 


0.02 


0.03 


0.03 


0.04 


0.05 


0.06 


0.07 


0.08 


0.10 


0.11 


0.12 


0.14 


0.15 


-3 


0.03 


0.04 


0.05 


0.06 


0.07 


0.09 


0.10 


0.12 


0.14 


0.16 


0.18 


0.20 


0.22 


-3+0.3 


0.04 


0.05 


0.06 


0.08 


0.10 


0.12 


0.14 


0.16 


0.19 


0.21 


0.24 


0.27 


0.29 


-3+0.7 


0.05 


0.07 


0.08 


0.10 


0.13 


0.16 


0.19 


0.22 


0.25 


0.28 


0.32 


0,36 


0.39 


-2 


0.06 


0.08 


0.10 


0.12 


0.15 


0.18 


0.22 


0.26 


0.29 


0.34 


0.38 


0.42 


0.47 


-2+0.3 


0.06 


0.08 


0.11 


0.14 


0.17 


0.20 


0.24 


0.28 


0.32 


0.37 


0.42 


0.47 


0.52 


-2+0.7 


0.06 


0.09 


0.11 


0.14 


0.17 


0.21 


0.25 


0.29 


0.34 


0.38 


0.43 


0.49 


0.54 


-1 


0.06 


0.08 


0.11 


0.14 


0.17 


0.20 


0.24 


0.28 


0.32 


0.37 


0.42 


0.47 


0.52 


-1+0.3 


0.05 


0.08 


0.10 


0.13 


0.16 


0.19 


0.22 


0.26 


0.30 


0.35 


0.40 


0.44 


0.49 


-1+0.7 


0.05 


0.06 


0.09 


0.11 


0.13 


0.16 


0.19 


0.23 


0.26 


0.30 


0.34 


0.38 


0.42 


0 


0.04 


0.06 


0.07 


0.09 


0.11 


0.14 


0.16 


0.19 


0.22 


0.25 


0.28 


0.32 


0.35 


0.3 


0.03 


0.05 


0.06 


0.07 


0.09 


0.11 


0.13 


0.15 


0.17 


0.19 


0.22 


0.24 


0.27 


0.7 


0.03 


0.03 


0.04 


0.05 


0.06 


0.07 


0.08 


0.09 


0.10 


0.11 


0.12 


0.14 


0.15 



Table 17. R for water- vapor window (660-1220 cm -1 ). 



Log it 


-80° 


-70° 


-60° 


-50° 


-40° 


-30° 


-20° 


-10° 


0° 


+10° 


4-20° 


4-30° 


4-40 0, 


-2 


0 


0 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


-2+0.3 


0.01 


0.01 


0.01 


0.02 


0.02 


0.02 


0.02 


0.03 


0.03 


0.03 


0.03 


0.04 


0.04 


-2+0.7 


0.03 


0.03 


0.04 


0.05 


0.05 


0.06 


0.07 


0.08 


0.09 


0.09 


0.10 


0.11 


0.12 


-1 


0.05 


0.07 


0.08 


0.09 


0.11 


0.12 


0.14 


0.15 


0.17 


0.18 


0.20 


0.22 


0,24 


-1+0.3 


0.09 


0.11 


0.14 


0.16 


0.19 


0.21 


0.24 


0.27 


0.30 


0.33 


0.36 


0.39 


0.42 


-1+0.7 


0.19 


0.23 


0.27 


0.32 


0.37 


0.42 


0.48 


0.54 


0.60 


0.66 


0.72 


0.78 


0.85 


0 


0.30 


0.37 


0.44 


0.52 


0.61 


0.70 


0.79 


0.89 


0.99 


1.09 


1.19 


1.29 


1.40 


0.3 


0.47 


0.57 


0.69 


0.82 


0.95 


1.10 


1.25 


1.40 


1.56 


1.72 


1.89 


2.05 


2.22 


0.7 


0.77 


0.95 


1.15 


1.37 


1.60 


1.84 


2.10 


2.36 


2.64 


2.91 


3.19 


3.48 


3.76 


+ 1 


1.00 


1.24 


1.51 


1.80 


2.11 


2.43 


2.77 


3.13 


3.49 


3.86 


4.24 


4.62 


5.00 


1.3 


1.18 


1.46 


1.78 


2.12 


2.49 


2.88 


3.29 


3.71 


4.14 


4.59 


5.04 


5.49 


5.95 


1.7 


1.26 


1.57 


1.91 


2.28 


2.68 


3.10 


3.54 


3.99 


4.46 


4.94 


5.43 


5.92 


6.42 
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Table 18. R for water vapor, total. 



Log u 


-80° 


-70° 


-60° 


-50° 


-40° 


-30° 


-20° 


-10° 


0° 


+10° 


+20° 


+30° 


+40° 


-5 


0.045 


0.053 


0.062 


0.073 


0.086 


0.103 


0.123 


0.148 


0.176 


0.208 


0.243 


0.282 


0.325 


-5+0.3 


0.066 


0.076 


0.088 


0.103 


0.120 


0.142 


0.168 


0.199 


0.234 


0.275 


0.320 


0.371 


0.426 


-5+0.7 


0.105 


0.119 


0.137 


0.158 


0.183 


0.214 


0.250 


0.292 


0.341 


0.397 


0.460 


0.530 


0.607 


-4 


0.146 


0.164 


0.188 


0.216 


0.248 


0.288 


0.334 


0.388 


0.451 


0.522 


0.602 


0.691 


0.790 


-4+0.3 


0.199 


0.223 


0.254 


0.291 


0.333 


0.384 


0.443 


0.512 


0.592 


0.682 


0.783 


0.896 


1.022 


-4+0.7 


0.291 


0.326 


0.370 


0.421 


0.480 


0.550 


0.632 


0.726 


0.834 


0.957 


1.094 


1.248 


1.419 


-3 


0.375 


0.421 


0.477 


0.541 


0.616 


0.704 


0.807 


0.925 


1.060 


1.213 


1.385 


1.577 


1.789 


-3+0.3 


0.470 


0.528 


0.597 


0.677 


0.771 


0.881 


1.009 


1.156 


1.323 


1.512 


1.725 


1.962 


2.223 


-3+0.7 


0.608 


0.684 


0.773 


0.877 


1.000 


1.144 


1.311 


1.502 


1.720 


1.967 


2.244 


2.552 


2.892 


-2 


0.715 


0.806 


0.912 


1.036 


1.183 


1.355 


1.553 


1.780 


2.040 


2.335 


2.665 


3.032 


3.438 


-2+0.3 


0.824 


0.931 


1.056 


1.202 


1.374 


1.575 


1.806 


2.072 


2.376 


2.721 


3.106 


3.535 


4.010 


-2+0.7 


0.975 


1.106 


1.259 


1.437 


1.645 


1.887 


2.166 


2.487 


2.853 


3.267 


3.729 


4.244 


4.815 


-1 


1.095 


1.246 


1.422 


1.627 


1.865 


2.141 


2.459 


2.824 


3.239 


3.707 


4.229 


4.812 


5.460 


- 1 +0.3 


1.224 


1.398 


1.600 


1.834 


2.105 


2.419 


2.780 


3.193 


3.661 


4.186 


4.772 


5.427 


6.155 


-1+0.7 


1.421 


1.633 


1.877 


2.159 


2.484 


2.857 


3.284 


3.769 


4.315 


4.926 


5.609 


6.369 


7.211 


0 


1.603 


1.852 


2.138 


2.468 


2.845 


3.275 


3.764 


4.316 


4.935 


5.626 


6.395 


7.248 


8.190 


0.3 


1.823 


2.120 


2.460 


2.850 


3.293 


3.794 


4.360 


4.995 


5.704 


6.491 


7.362 


8.324 


9.382 


0.7 


2.177 


2.554 


2.985 


3.474 


4.027 


4.649 


5.345 


6.120 


6.977 


7.921 


8.957 


10.091 


11.329 


+ 1 


2.431 


2.869 


3.368 


3.933 


4.569 


5.282 


6.077 


6.959 


7.930 


8.993 


10.153 


11.416 


12.787 


1.3 


2.610 


3.092 


3.640 


4.261 


4.960 


5.743 


6.613 


7.576 


8.634 


9.790 


11.047 


12.410 


13.884 


1.7 


2.686 


3.182 


3.751 


4.397 


5.127 


5.946 


6.858 


7.866 


8.973 


10.182 


11.497 


12.922 


14.461 



Table 19. —dR/d log u for water vapor, total. 



Log u 


-80° 


-70° 


-60° 


-50° 


-40° 


-30° 


-20° 


-10° 


0° 


+10° 


+20° 


+30° 


+40° 


-5+0.3 


0.08 


0.09 


0.10 


0.11 


0.13 


0.15 


0.17 


0.19 


0.22 


0.26 


0.29 


0.33 


0.38 


-5+0.7 


0.12 


0.13 


0.15 


0.17 


0.19 


0.22 


0.25 


0.28 


0.32 


0.37 


0.42 


0.47 


0.54 


-4 


0.16 


0.17 


0.19 


0.22 


0.25 


0.28 


0.32 


0.37 


0.42 


0.48 


0.54 


0.61 


0.69 


-4+0.3 


0.20 


0.22 


0.25 


0.28 


0.32 


0.36 


0.41 


0.47 


0.53 


0.60 


0.68 


0.77 


0.86 


-4+0.7 


0.26 


0.29 


0.33 


0.37 


0.42 


0.47 


0.54 


0.61 


0.69 


0.78 


0.89 


1.00 


1.12 


-3 


0.30 


0.34 


0.38 


0.43 


0.49 


0.55 


0.63 


0.72 


0.82 


0.93 


1.06 


1.19 


1.34 


-3+0.3 


0.33 


0.37 


0.42 


0.48 


0.55 


0.62 


0.71 


0.81 


0.93 


1.06 


1.21 


1.37 


1.55 


-3+0.7 


0.35 


0.40 


0.45 


0.52 


0.60 


0.69 


0.79 


0.90 


1.04 


1.19 


1.36 


1.55 


1.76 


-2 


0.36 


0.41 


0.47 


0.54 


0.62 


0.72 


0.83 


0.95 


1.10 


1.26 


1.44 


1.64 


1.86 


-2+0.3 


0.37 


0.43 


0.49 


0.57 


0.65 


0.75 


0.87 


1.00 


1.15 


1.32 


1.51 


1.72 


1.95 


-2+0.7 


0.39 


0.46 


0.53 


0.61 


0.70 


0.81 


0.94 


1.08 


1.24 


1.42 


1.61 


1.83 


2.08 


-1 


0.41 


0.49 


0.57 


0.66 


0.76 


0.88 


1.02 


1.17 


1.34 


1.53 


1.73 


1.96 


2.22 


-1+0.3 


0.45 


0.54 


0.63 


0.73 


0.85 


0.99 


1.14 


1.31 


1.49 


1.69 


1.91 


2.16 


2.44 


-1+0.7 


0.55 


0.66 


0.78 


0.92 


1.08 


1.25 


1.43 


1.63 


1.85 


2.10 


2.37 


2.65 


2.95 


0 


0.67 


0.81 


0.97 


1.15 


1.35 


1.56 


1.79 


2.04 


2.30 


2.59 


2.91 


3.25 


3.62 


0.3 


0.81 


0.99 


1.19 


1.41 


1.65 


1.91 


2.19 


2.50 


2.83 


3.19 


3.58 


3.99 


4.43 


0.7 


0.90 


1.10 


1.33 


1.59 


1.87 


2.18 


2.52 


2.88 


3.26 


3.66 


4.08 


4.52 


4.98 


+ 1 


0.74 


0.92 


1.13 


1.36 


1.61 


1.89 


2.19 


2.51 


2.85 


3.21 


3.59 


4.00 


4.43 


1.3 


0.40 


0.51 


0.63 


0.76 


0.90 


1.06 


1.24 


1.44 


1.66 


1.90 


2.16 


2.44 


2.74 


1.7 


0.06 


0.06 


0.06 


0.07 


0.09 


0.12 


0.16 


0.21 


0.27 


0.34 


0.42 


0.51 


0.61 



Table 20. Chart areas, fRdT at -80° (u = const.). 



Log M 


co 2 


Log u 


o 3 


H2O total 


6.3 m band 


Rotat. band 


Window 


-4 


0.18 


-5 




3.72 
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